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By Doris Cohen 


SUMMARY 

The method oj superposition of linearized conical flows has 
been applied to the calculation of the aerodynamic properties ^ in 
supersonic flighty of thin flat ^ swept-hack wings at an angle of 
attack. The wings are assumed to have rectilinear plan forms ^ 
with tips parallel to the stream, and to taper in the conventional 
sense. The inpestigation covers the moderately supersonic speed 
range where the Mach lines from the leading-edge apex lie ahead 
of the wing. The trailing edge may lie ahead of or behind the 
Mach lines from its apex. The case in which the Mach cone 
from one tip intersects the other tip is not treated. 

Formulas are obtained for the load distribution, the total lift, 
and the drag due to lift. For the cases in which the trailing edge 
is outside the Mach cone from its apex {supersonic trailing edge), 
the formulas are complete. For the wing with both leading and 
trailing edges behind their respective Mach lines, a degree of 
approximation is necessary. It has been found possible to give 
practical formulas which permit the total lift and drag to be 
calculated to within 2 or S percent of the accurate linearized- 
theory value. The local lift can be determined accurately over 
most of the wing, but the trailing-edge-tip region is treated only 
approximately. 

Charts of some of the functions derived are included to facili- 
tate computing, and several examples are worked out in outline. 

INTRODUCTION 

It is customary, in supersonic wing theory, to describe 
any straight segment of the boundary of a wing plan form as 
supersonic or subsonic accordingly as the segment lies out- 
side or is contained within its foremost Mach cone; that is, 
as the component of the flight velocity normal to the edge is 
greater than or less than the speed of sound. These two 
circumstances result in fundamentally different types of flow 
over the surface. It is apparent that the real reference is 
not to a property of the wing plan form, but to a combination 
of plan-form geometry and the velocity of the wing relative 
to the speed of sound. Thus (see fig. 1) every swept-back 
wing, on entering the supersonic regime, has subsonic leading 
and, in most cases, subsonic trailing edges. At a higher Mach 
number, the same plan form may have subsonic leading edges 
and supersonic trailing edges. Finally, if the Mach number is 
increased suflSciently, both leading and trailing edges will 
become supersonic. 


Interference effects also depend on the flight Ma<jh number, 
since the extent of the various disturbance fields is determined 
by the angle between the Mach lines. Thus, no single 
concise formula or method of treatment has as yet been 
developed to predict, even approximately, the aerodynamic 
characteristics of an arbitrary wing plan form through the 
supersonic speed range. 

The present report is concerned with the loading, lift, and 
drag, according to linearized theory, of thin, flat, swept- 
back wings with rectilinear boundaries and conventional 
taper. Various methods are available for the calculation 
of these properties when the leading edge is supersonic. 
Of these, the method of reference 1 is perhaps the most 
convenient. Formulas obtained by this method for the 
loading and lift-curve slope of wings with supersonic lead- 
ing and trailing edges are presented in reference 2. In the 
following, therefore, the emphasis will be on the solution of 
the problems arising from the interaction of the flow fields 
in the presence of subsonic leading edges (figs. 1 (b), (c), 



Figure l.— A typical tapered swept-back wing at six supersonic Mach numbers, showing 
the Mach lines from the leading- and trailing-edge apexes and from the tips. 


and (d)). The case (fig. 1 (a)) in which the Mach number 
and aspect ratio are so low that interaction takes place 
between the tip flow fields will not be treated. An approxi- 
mate solution to this problem may be found in reference 3. 

When a wing with a subsonic leading edge is to be studied, 
considerable simplification of the problem may be achieved 
by making use of the solutions, available in reference 4 

1 


2 
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and other sources, for the inhiiite triangular wing.^ From 
these solutions the aerodynamic characteristics of a variet}^ 
of swept-back plan forms can be calculated by the use of 
the superposition principle of linearized theory to cancel 
any lift beyond the specified wing boundaries. Two methods 
of cancellation have been developed: one, presented in refer- 
ence 5, uses supersonic doublets and is general enough to 
apply to curved boundaries; the other, originally due to 
Busemann (reference 6), cancels by means of the super- 
position of conical flow fields. In the present report the 
conical-flow method is used, since it appears to offer some 
advantages for the straight-sided plan forms under 
consideration, particularly in determining the integrated lift. 

The material presented in this report is largely drawn 
from references 7, 8, and 9, with some simplifications sug- 
gested by practical exjierience. In particular, the formulas 
for the total lift hav(i been reworked to substitute, with 
no increase in computational labor, a combined “primary'^ 
and ^^secondary'^ correction for each of the ^‘primary^^ cor- 
rections in reference 7. Also, the formulas containing elliptic 
integrals have been rewritten to take full advantage of 
available tables. As in the preceding papers, the final for- 
mulas will be derived for unyawed wings with tips parallel 
to the stream, but the application of the general method 
and the basic solutions to other plan forms and problems 
will be apparent. Some numerical examples will be included 
in order to show the magnitude of the effects discussed and 
to summarize the method. A table summarizing the 
formulas is also included. 

I— METHOD OF THE SUPERPOSITION OF CONICAL 

FLOWS 



A conical flow field is one in which the velocity components 
u, and w in the stream, cross-stream and vertical directions, 
respectively, are constant in magnitude along any ray from 
the foremost point, or apex, of the field. Such flows are 
found as solutions of the linearized potential equation for 
supersonic flow. A detailed discussion of their derivation 
and use is contained iu reference 4. In the cancellation-of- 
lift procedure, only solutions of the so-called “mixed^’ type 
described in section V of reference 4 are required, except for 
the basic solution (for the infinite triangular wing) which is 
itself of conical form. 

SYSTEM OF NOTATION FOR CONICAL FLOWS 

The Cartesian coordinate system is placed so that the 
origin coincides with the projection of the leading-edge apex 
on the horizontal plane , the positive x axis extending down- 
stream from the origin and the y axis extending perpendicular 
to the X axis in the horizontal plane. (See fig. 2.) For the 
conical flow fields, it is further convenient to define a variable 
to designate a particular ray in the xy plane, since the flow 
velocities are constant along such a ray. If the apex of the 
field is specified, then tlie ray is most readily described by its 
slope, measured from, the downstream direction. The 
conical solutions of the. supersonic flow equation are, how- 

* The present report covers in detail only unyawed wings. However, yawed wings may 
be treated similarly, starting with the yawed triangular-wing solutions. This problem 
is the subject of a paper, NACA 1'N 2262, 1950, by Lampert, prepared concurrently with 
the present report. 


ever, functions of the ratio of the slope of the ray to the 

slope i of the Mach lines, where is and M is the 

P 

free-stream Mach number. For the triangular- wing flow 
with its origin at the apex of the wing, therefore, the conical 
variable will be chosen as 


a=^ 


y 

X 


( 1 ) 


At the Mach lines from the leading-edge apex, a equals ±1. 


The ray from 0, the wing apex, making the angle tan“^ — 


with the stream will hereinafter be referred to as the ray a, 
and the subscript a will indicate association with a constant- 
load sector (to be introduced later) of which such a ray is 
one of the boundaries. 

For each of the conical fields to be superposed at the* edges 
of the wing plan form, a new coordinate system is set up with 
its origin at the apex of the field. In conformity with the 
notation of reference 4, the conical variable relative to the 
displaced origin is called <, with subscripts to denote the 
location of the origin. Thus, if Xa^ya is the point of intersec- 
tion of the ray a with the plan-form boundary and is to serve 
as the apex of a canceling conical field. 




y-Va 

X — Xa 


( 2 ) 
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is the ratio of the slope of the ray ta of that field to the slope 
of the Mach liaes. 

If the ratio of the slope of the leading edge to the slope of 
the Mach lines is 

m=/3 cot A (3) 


(4) From equation (5), u must equal zero in the portion 
of the horizontal plane not covered by conditions (2) or (3). 

(5) In supersonic flow there exists the additional condition 
that all the velocities must go to zero on the Mach cone from 
the apex of the field. 


where A is the angle of sweepback, then at the leading edge 
a=m, and a ray from the leading-edge tip is designated by 
tm- If s is the wing semispan, the leading-edge tip has the 
3 s 

coordinates — ; s and any point x^y has the conical coordinate 


t =3~ y~ ^ 
^ 3s 

X 

m 


(4) 


in the field with apex at — ; s. 

m 

I 

Other symbols referring to angular locations will be defined 
in the same way as needed. A summary of the symbols 
will be found in appendix A. 

BOUNDARY CONDITIONS FOR CANCELLATION OF LIFT 

The general problem of deriving the flow over a wing of 
finite dimensions from the known flow over an infinite wing 
is the problem of determining the induction effects due to 
the edges. These effects may be thought of as associated 
with the cancellation of the lifting pressure at the boundaries 
of the finite wing. In the linearized lifting-surface theory, 
they may be evaluated by the superposition of flow fields 
with negative lifting pressure over the portion of the infinite 
wing outside the boundaries of the finite plan form, provided 
the other boundary conditions are not disturbed. In the 
case of a flat wing at an angle of attack, the latter provision 
means that the canceling field must (1) induce no down wash 
within the boundaries of the finite wing and (2) introduce no 
new lifting pressure outside those boundaries. 

In accordance with thin-airfoil theory, the boundary con- 
ditions will be satisfied in the horizontal jflane rather than on 
the surface of the wing. Also, by thin-airfoil theory, the 
conditions on the lifting pressure are converted to conditions 
on the velocity field through the relation 

(?l.+o 

In the simplest case, the lift to be canceled will be dis- 
tributed uniformly over a semi-infinite region bounded by 
two straight lines. The boundary conditions of the problem 
may then be said to be conical with respect to the intersection 
of the two lines, which become “rays’' of the canceling conical 
field. The boundary conditions on the canceling velocity 
field in this case may be summarized as follows : 

(1) The stream wise velocity u must approach values equal 
in magnitude and opposite in direction on the upper and lower 
surfaces of the horizontal plane. 

(2) In the horizontal plane, u must be constant over the 
infinite sector in which lift is to be canceled. 

(3) The vertical velocity w must be zero in the portion 
of the 2=0 plane occupied by the projection of the finite 
wing. 


CANCELLATION OF NONUNIFORM LIFT 

The foregoing are the general conditions for a uniformly 
loaded canceling flow field. Under the proper conditions, a 
nonuniform distribution of lift may be canceled by the super- 
position of a number of such fields. This procedure is best 
explained by a concrete example. 

Consider the problem of a swept-back wing flying at a 
high Mach number such that, as in figure 1 (e), the Mach 
lines from the leading-edge apex intersect the tips of the 
wing. The method of deriving the swept-back wing from 
an infinite triangular wing in that case is indicated in figure 
3. It may be noted at the start that, according to linear 
theory, the lift behind the supersonic trailing edge may be 
canceled in any way without affecting the velocities on the 
wing. Thus it remains only to consider the effect of can- 
celing the lift outboard of the tips. 



Figure 3. — Method of cancellation of lift beyond the tip when the leading-edge Mach 
line intersects the side edge of the wing. 


An infinite triangular wing with supersonic leading edges 
has a load distribution which is constant over the portions 
of the wing between the leading edge and the Mach lines 
from the leading-edge apex (see fig. 4). This constant load 
may be canceled outboard of each of the tips of the swept- 
back wing by a single negatively loaded triangle of infinite 
extent, one side coinciding with the side edge of the wing 
and a second side coinciding with the extension of the leading 
edge. However, the area to be removed (region BAG, fig. 
3) includes also a region over which the pressure varies, and 
is conical with respect to 0. Since the boundaries of the 
region are conical with respect to A, no one conical solution 
can satisfy the requirements of the problem. The problem 
is brought within the limitations of the conical solutions by 
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region of uniform loading 

Figure 4.— Lift distribution on a triangular wing with supersonic leading edges. 


considering the lift to be made up of an infinite number of 
constantly loaded, ovei'lapping sectors of infihite extent. 
(See fig. 3.) These sectors are bounded on one side by the 
wing tip ; the second sidfi is the extension of a ray from apex 
0 of the wing. Between the leading edge {a=m) and the 
leading-edge Mach line (a=l), no division of the field is 
necessar}^ since the lift density is constant in that region. 


7Yt 

If a sector with apex at .k and angle tan“^— is used to cancel 


this uniform lift, then the remaining superposed fields must 
be used where a<^l (see fig. 4) to restore the difference 
between that lift and the loading on the triangular wing. 

If Ui is the stream wise component of the perturbation 
velocity corresponding to the constant loading ahead of the 
leading-edge Mach lines, and u^{a) is the same velocity in 
the region between the Mach lines, then the magnitude of 
the u component of the velocity in the initial canceling 
sector will be —Ui, and on the remaining sectors (see fig. 3) 
minus the increment in corresponding to an increment 

in a, or da. (Note that this last quantity is positive, as 

required). To determine the total effect of canceling the 
loading outboard of th<5 tip, the velocities induced by the 
latter infinitesimally loaded elements are integrated and 
added to the negative effect of the initial constant-load 
sector. 


II— LOADING ON WING WITH SUBSONIC LEADING 

EDGE 


LOAD DISTRIBUTION OVER TRIANGULAR WING 


In the notation of this paper, the velocity distribution 
over a flat lifting trianglB with leading edge behind the Mach 
lines may be written 




muQ 

■ylrn^—a^ 


(6) 


where 


mV a 

(m) 


(7) 


is the (constant) velocity along te center line a = 0. In the 
expression for Uoj E'im) is th e comp lete elliptic integral of 
the second kind, of modulus The load distribution 

is obtained from the velocity distribution by equation (5). 

SWEPT-BACK WING WITH SUPERSONIC TRAILING EDGE 
(TIP CORRECTION) * * 

If the problem is now to find the loading on a swept-back 
wing with subsonic leading edges, but supersonic trailing 
edges, only the tip effects will modify the triangular-wing 
distribution. The calculation of the tip effect on a wing 
with subsonic leading edge (m<l) is somewhat complicated 
by the fact that the pressure becomes infinite at the leading 
edge, but otherwise follows the procedure outlined in the 
preceding section. 

It will first be necessary to present the expression for the 
previously described conical field with uniformly loaded sec- 
tor to be used as the element in canceling the lift outboard 
of the tip. 

ELEMENTARY SOLUTION FOR A STREAMWISE TIP 

If s is the semispan of the wing, the apex of any element a 
(see the section on Notation) is at 



and, from equation (2), 



a 


Then, if Ua is the constant perturbation-velocity component 
to be canceled over the region between the tip and the exten- 
sion of the ray a, the previously listed boundary conditions 
for each of the required canceling fields may be written as 
follows (see fig. 5): 

(1) and (2) When 0<<a<c^, u=^±Ua (constant for the field) 

(3) When ^^<0, w=0 

(4) When u=0 

(5) When u=v=w=(). 

The solution of the supersonic flow equation satisfying the 
above boundary conditions has been derived in reference 4.^ 
In the xy plane, the stream wise component of the velocity is 


, Ua _i a-\-ta~\~^ata 
u=±r.^. — cos — 7 

7T ta a 


( 10 ) 


The signs refer to the upper and lower surfaces, respectively. 

In figure' 5, the essential features of the solution are 
shown. At the top is a detail view of the wing side edge and 
shows the boundary conditions. In the center is a typical 
plot of the argument of the inverse cosine in equation (10), 
against Where this quantity is less than —1 (i. e., 
0<fo<a), the real part of the inverse cosine is ir. Where 
the argument is greater than +1 (<o!i>® <a<C~l)> the 

2 Approximate formulas, valid when m is close to 1, have been presented for this case in 
reference 10. 

* The corresponding solutions for raked-in or raked-out tips may also be found in reference 4, 
or deduced from later sections in the present report. 
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ia 

Figure 5.— Elementary solution for canceling lift at the tip. 

real part of the inverse cosine is zero. On the wing 
(— the argument goes from +1 to —1 and the 
inverse cosine is real. Thus in canceling, or subtracting, 
the velocity Ua between <a=0 and ta=a^ the increment in 
velocity 

'I/ n -\—i -i-Oni 

( 11 ) 


^ (a:, y, a)= ^ cos ^ ^ 

7T tn a 


is induced on the wing upper surface. 

TIP-INDUCED CORRECTION TO THE LOADING 

Following the procedure outlined in Part I, we proceed to 
determine the effect of canceling the lift outboard of the 
wing tip. Since the value of Ua for the initial canceling 
field —u^{m) and the value for the first incremental field 

da are both infinite when the leading edge is subsonic, 

it is first necessary to write the induced velocity at a point 
x,y as 


(Au) 


tip 


976745-5^- 


_ lim r cos“^ a+ta-\-2atg 

a->m L 7T - 4 ^ “ 


ta—a 


i p ^ cos-1 da] (1 2) 

wjao da tg—a J ^ 


where the limit a^ is the value of a corresponding to the 
rearmost sector including the point x,y in its Mach cone. 
The value of a^ is found by setting tg (equation (9)) equal 
to —1. Thus, for the tip correction, 


O-o — 




^+&{y-s) 


(13) 


This parameter will be additionally useful as the value of a 
at which the velocity" correction given by equation (11) 
goes to zero and its derivative has a singularity. 

Before performing the integration of equation (12), tg 
must be replaced by its expression in terms of x, y, and a. 
Then integration of the second term by parts results in a 
term which, at the upper limit, exactly cancels the first 
term, and at the other limit is zero, leaving, after substitution 
for Ui,, 

— m(a^ + /3y)-Uo f"* Vao(.9 — y) da 

J«0 (ax-^y)^(^2_a2) (1 +a) {a-a^) 

(14) 

This integral is finite and can be evaluated in terms of 
elliptic integrals as follows: 


where 


Ao=KoE{4^, k')-{Ko-Eo)F{^, r) 


(16) 


and Kq and Eq are 2/tt times the complete elliptic integrals 
K and E of modulus 

, ^ / (m — ao)(l— m) 

V 2m(ao+l) 

In equation (16), F and E are the incomplete 

integrals with the complementary modulus k' = ^|\—k‘^ and 
argument 




/3s(ao+m) 

The functions Kq, Eq and Aq are tabulated in reference ! 1 ^ or 
ma}^ be computed from the tables of reference 12. A plot 
of Ao is given in figure 6. 

Value at the side edge. — At the tip, y is equal to s and the 
first term in equation (15) vanishes. In the second term, ^ 
becomes tt/2 and E {\k,E) and F (^,/r') reduce to the complete 
integrals E'=E{k') and K'=K{k')y respectively. Then, 
since, by Legendre’s relation, 

K'E-K'K+KE'=wl2 

Ao reduces to 1. The induced velocity correction is seen to 
be exactly equal to —Uaj bringing the lift to zero at the 
wing tip. 

Drop in lift across tip Mach line. — An interesting effect 
shows itself at the other limit of the tip region, that is, at the 
Mach line from the tip of the leading edge. Along this line 
only the influence of the leading-edge pressure is felt, so that 


* The quantity Kc is called Fo in reference 11. 
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Figure 6— The elliptic function Ao(^,k)=KoEi4^,k')—(Ko—ED)Fi^X). 
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ao=m. Tlien k=0, k'=l, K=ttI 2, E{yp,k') reduces to 
sin ^ = -y/ finally 


(A«),(p(<„ = — l)=Aw* = 


— Wog.9 

A 2(z + /3i/)(mx — /3«) 


(17a) 


or, since along the lip Mach line /3(s— ?/) = 


— y 

m 


Au* = 


—UQ\m^s 
-\'2(1 + m){mx — ps) 


(17b) 


Tills result indicates a finite drop in pressure across the 
Mach line from the tip, an effect which is associated with the 
cancellation of infinite pressure at the leading edge and con- 
sequently does not appear as long as the leading edge is 
ahead of the Mach lines. The ratio of the drop in lift across 
the tip Mach line to the uncorrected lift can be written 


Au* _ 1 (1 -fa) (m-\-a) 

V 2m(l+m) 


(18) 


This ratio is plotted against a/m in figure 7 and shows the 
percentage loss of lift at the tip to be very large. In fact, 
for any but the lowest-aspect-ratio wings, the lift remaining 
in that region is almost negligible. This effect, which should 
be of considerable practical interest, was first indicated in 
the results of reference 13 for the limiting case of m=0. 



SWEPT-BACK WING WITH SUBSONIC TRAILING EDGE 

The tip-effect correction just derived applies equally to 
wings with supersonic or subsonic trailing edges. The effect 
of a subsonic trailing edge is calculated separately, and is 
primarily due to canceling the triangular-wing leading in 
the wake region. If, however, the triangular-wing loading 


has been modified by the introduction of side edges, then this 
modification must also be taken into account when canceling 
the lift behind the trailing edge. In the conical-flow method, 
the various component flow fields must be canceled individu- 
ally. The sections immediately following will discuss the can- 
cellation of the triangular-wing loading ; cancellation of the tip- 
induced components of velocit}^ will be considered under the 
heading “Secondary Corrections.’^ 

PRIMARY TRAILING-EDGE CORRECTIONS 

Procedure for canceling lift in the wake region. — The 
basic procedure is again to consider the load to be canceled 
to be built Up by the superposition of uniformly loaded 
sectors, bounded on one side (see fig. 8) by the rays a, and 
on the other by the trailing edge of the wing. It is con- 
venient at this point to introduce the parameter 

^t=/3Xcot (angle of sweep of trailing edge) 



Figure 8.— Oblique constant-lift element (shaded) for cancellation of lift at subsonic 
trailing edge, and induced velocity distribution. 


The boundary conditions to be satisfied by the u compo- 
nent of the elementary canceling velocity field are indicated 
for the right span in figure 8; each field must have constant 
velocity Ua when a<ta<rrit and zero stream wise velocity over 
the wake region, The concomitant vertical 

velocity must be zero on the wing surface. However, when 
a is small, the region — I <^a<a will ipclude a portion of the 
left-hand wing panel. Since in this region the u component 
of velocity has already been specified, the vertical velocity 
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will not, in general, be zero. Nor is it possible to modify 
the field to satisfy the boundary condition on the far wing, 
since the area involved is not conical with respect to the apex 
of the field. 

The error involved in i:he foregoing procedure is minimized 
by the use of a symmetrical flow field to cancel the initial 
load Uq at a=0, where a single conical field can be made to 
satisfy the boundary conditions exactly on both wing panels. 
This flow field (see fig. 9) would have its origin at the apex 
Co,0 of the trailing edge, and the constant-load region would 
extend over the entire wake region. Between the trailing 
edge of the wing and the Mach lines from Co,0 the induced 
downwash would be zero in the plane of the wing, w^hile the 
pressure would vary as required to satisfy the fundamental 
flow equations. 

In figure 9, a typical curve of is shown, from which it 
can be seen that the load to be canceled is very nearly 
constant over a considerable fraction of the wake region. 
Cancellation of the velocity by the symmetrical field will 
consequently leave only a small variation in u to be canceled 
by the oblique fields d(»scribed earlier in the section. The 
resulting violation of the; flat-plate condition may be expected 
to be small, ^ and wdll take place only over a small region 
near the tip of the trailing edge. 



Figure 9.— Symmetrical field for cancellation of wo at subsonic trailing edge. 


Symmetrical solution.— For the symmetrical solution we 
define the conical variable 


— 


X—Co 


(19) 


which is zero along the center line of the wing and equals ±mt 
at either trailing edge. Then the boundary conditions to be 
satisfied in the xy plane may be summarized as follows: 


--mt<to< + 7nt 

rritK'M ^ 1 


u= ±Uo 
w = 0 


The required solution is given in reference 14. The u com- 
ponent in the xy plane is 

where K'{m^ is the complete elliptic integral of the first kind 
of modulus \\—m} and is the corresponding 

incomplete integral of argument 

/ 1 _/ 2 

i • _ 1 / J ‘0 

0 = sin \ 2 

VI —m/ 

The form of the induced velocity on the wing (see fig. 9) is 
very similar to the inverse cosine curves of the tip solutions. 

On the wing, <t> is real and the symbols r. p. may be omitted. 
The velocity induced on the upper surface by cancellation of 
Uo behind the trailing edge is therefore 


{Au)o= - 


Uq 

K'{m,) 


F(</>,^ I— mr^) 


( 20 ) 


* Calculations made to check this statement have shown the induced downwash angles to 
be less than 0.5 percent of the angle of attack, even in the most unfavorable circumstances. 


Oblique solutions for the wake region. — The symbol 
ta will be used as before to indicate a ray of the flow field 
with apex at the point of intersection of the ra\’ a with 

the wing boundary — in this case the trailing edge. Along 
the trailing edge, 


rrit , 

Va=-^{Xa 


-Co) 


Since a = 0 (ValXa), we may solve for Xa and ?/« as functions 
of a and the constants m, and Cq: 


Then 


rritCo 

Xfi 

(21) 

rrit — a 

. rritCoa 

(22) 

^ fiy(m, — a) — m,Coa 
“ x(mt—d) — 7ntCo 

(23) 


The boundary conditions to be satisfied by the elementary 
solution are (for a>0) 


a<ta<mt 
mtKta< + l 

-l<ta<a 


U= ±Ua 
W=0 
U = 0 
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solution satisfying these conditions may be obtained 
from the tip solutions by an oblique transformation. (See 
reference 4.) In the xy plane, the resulting expression for 
the u component is 

Then the velocity induced at any point on the upper 
surface of the wing by the cancellation of the infinitesimal 
increment of perturbation velocity Ua over the sector bounded 
by the ray a and the trailing edge is 

^ = cos- n_-«) (I 

(fa 7T {\—m,)(ta—a) 

(24) 

Correction of loading near the trailing edge. — To deter- 
mine the lift at any point Xjy near the trailing edge of the 
wing, it is first necessary to determine the most rearward 
canceling sector Gq that will influence that point. Setting ta 
(equation (23)) equal to 1, we solve for 


^ _ x-py-Co 

Gq — TTl t a 

x — Py — mtCo 


(25) 


Then the total correction to the triangular-wing velocity 
obtained as a result of canceling that velocity behind the 
trailing edge is 

r^o rlAv 

(Au)r.£. (x, 2 /) = (Am)o+J^ da (26a) 

The integral in the foregoing expression has been evaluated 
in terms of an incomplete elliptic integral of the third kind, 
which may be computed with the aid of the tables of refer- 
ences 1 1 and 15. Because it will be necessary to define several 
new functions it was thought better to present the results in 
an appendix (appendix B). For practical use, graphical or 
numerical integration may be preferred, in which case a 
convenient form is obtained by rewriting Ua as {du^jda) da^ 
or duA, in equation (24). Thus equation (26a) becomes 


{Au)t.e. {x,y) = {Au)o— 


1 

- cos"' 

7T Jwo 


(1— a)(fg— mj — (m< — g)(l— g 

(1— mO(f«— a) 


duA (26b) 


where ta and Ua must be evaluated for selected values of a 
between zero and Gq. The integrand, of course, goes to zero 
at Ua(go)- At points along the leading edge (in cases in which 
the leading edge extends into the zone of influence of the 
trailing edge), the integral takes on a somewhat simpler form, 
with the result that the entire trailing-edge correction at 
such points can be written 


(Au)r..\x, «oj 

1 , 

LOS . V I 

7T 




(26c) 


where the first term inside the braces is and, in the last 

Uo 

term, 

/ • _i / 2^10 

\p=sm \ — j — 

\ m + Go 

and 

, ^ I jl—m) (m + Go) 

V 2m(l— Uo) 

SECONDARY CORRECTIONS 

The term “secondary corrections^^ is used here to designate 
the effect of canceling the lift introduced outside the bound- 
aries of the wing in the process of canceling the original tri- 
angular-wing loading beyond the tips and behind the trailing 
edge. As previously mentioned, cancellation of lift at the 
tip introduces new (negative) components of lift to be can- 
celed at the trailing edge. The original cancellation of lift 
behind the trailing edge, on the other hand, will introduce 
negative incremental pressures outboard of the tip and, under 
certain circumstances (see figs. 1 (a) and (b)), ahead of the 
leading edge. The distribution of lift to be canceled in each 
case is no longer part of a single conical field, but is composed 
of an infinite number of superposed conical fields originating 
at various points along the trailing edge or tip. In order to 
cancel these pressures accurately, it would be necessary to 
set up, for each of the original canceling elements, an infinity 
of positively loaded elements at the opposite boundary. 
Thus, each secondary correction would require a double 
integration for each point, and would obviously be quite 
tedious. The procedure is described in detail in references 
7 and 8. The more recent work of Mirels (reference 5) offers 
an alternative method which, while no less tedious at the com- 
putational state, is somewhat easier to set up for computing. 
Nevertheless, the exact calculation of the secondary correc- 
tions, and of the succeeding corrections arising as the second- 
ary corrections are in turn canceled at the opposite edges, 
appears feasible only with the aid of high-speed computing 
machinery. 

These corrections may be thought of as a converging 
series, since in each case (ex(5ept in the neighborhood of the 
leading edge) the induced effect is smaller than the canceled 
lift. Over most of the wing, the secondary correction is of 
the same order of magnitude as the tolerable error. Formulas 
for obtaining a major part of the secondary corrections can be 
given rather simply and should suffice to give results of 
practical accuracy in problems (fig. 1(c)) not involving lead- 
ing-edge corrections. Problems of the type shown in figure 
1 (b) will be discussed in a later section. 

Secondary corrections at the trailing edge. — The pressure 
differences induced by the tip are in the main due to can- 
cellation of the infinite pressure at the leading edge. It 
should therefore be permissible, for the secondary corrections, 
to approximate the tip-correction field by a single conical 
field from the leading-edge tip. The lift associated with this 
field may then be canceled behind the wing (see fig. 10) by a 
single infinity of superposed fields, as was the original triangu- 
lar-wing loading. If the values of (Au)tip calculated for 
points xif.yi, along the trailing edge are assumed to apply all 
along the corresponding ray^s tm(xt,yb) from the tip, then the 
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Figure 10 —Sketch for approximate cancellation of extraneous lift introduced behind 
the hjading edge by the tip correction. 


lifting pressure will be exactly canceled along the trailing edge 
and the remaining variation of pressure in the wake will have 
very little effect on the flow over the wing. 

The cancellation fields are of the previously used oblique 
type, with a replaced by 


tm 


0(yb—s) 



(27) 


Let the particular point at which the line tm= — l intersects 
the trailing edge be designated by and other symbols 

referring to that point be similarly starred. Then the 
velocity induced at any point x^y on the wing by removal of 
{Au)tip along the trailing edge will be (from equation (24)) 


—Au* — 


I rtm(^o,vo) d{Au)tip 
TtJ-I dtm 


COS 




dtm 

(28) 


where Au* is given by equation (17) ,t* and ti, are calculated by 


A* &(y-y*) 


(29) 


and 


U 


P( y-yb) 

x — x^ 


( 30 ) 


respectively, and Xo,yo is the point of intersection of the Mach 
forecone from /, y with the trailing edge. 

The derivative (Au)tip would have to be determined 
atm 

numerically or graphically from a plot of the calculated 
values of (Am) u„ against t„. In order to avoid this procedure, 
it is preferable to rewrite expression (28) as 

—Am* 2(t* — m,) — {m, + l){l—t*) 

— (i_w, )(<*+!) 

1 f(Au),.p(xo.!/o) _ 

— I COS jz \ / . U\^y')tip 

TtJau* (1 — T^t)(tb — tm) 

(31) 


and integrate by plotting the inverse cosine function against 

(AV/) tip. 

As long as the aspect ratio of the wing is greater than 1//3 
(a condition already imposed by the exclusion of the problem 
ghown in fig. 1 (a)), Au* will be more than half (Au)tip at 
any other point on the trailing edge. Since, moreover, the 
integral term in equation (31) has zero slope at the Mach 
cone (^*=1), while the first term starts with infinite slope, 
it IS apparent that the secondary correction may be simplified 
still further by omitting the calculation of the integral. 
For points near the trailing edge, the loading can usually be 
faired to zero with sufficient accuracy. 

Secondary correction at the tip. — A similar method of 
approximating the secondary correction at the tip cannot be 
formulated with equal confidence. Since, however, over 
most of the wing the s^unmetrical correction (Ati)o contrib- 
utes the larger part of the total subsonic-trailing-edge 
effect, it will again be assumed that the entire effect con- 
stitutes a single conical field, with its apex at Cq, 0. The u 
velocity along each ray U will have the value (Au)t.e. (^by s) 
of the trailing-edge correction at the intersection s of the 
ray with the tip. The canceling fields will have the same 
form as those (equation (10)) used in deriving the primary 
tip correction, and the total approximate correction to the u 
velocity will be 


I riAu] 


(Au) 7*^ (Xq, 8 ) 


COS — diAujT. E. 

tb — ^0 


in which 


tb — 


0(y—s) 

X — Xb 


(32) 

(33) 


Xo is the value of Xb which makes — 1, and (Au)t.e. is 
calculated for x=Xb, y=s by equation (26b). 

NUMERICAL EXAMPLE 


Before proceeding to consider the problem of interaction 
between the leading and trailing edges, which introduces 
some radically different effects, the results so far obtained 
will be illustrated by a numerical example. The loading 
over an untapered wing, with cot A = 0.6 and reduced 
aspect ratio /3^=1.92, has been calculated at four spanwise 
stations: 25-, 50-, 75-, and 95-percent semispan. The wing 
plan form and section lift distributions are shown in figure 11. 


x — x 


0^0 
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(a) Section A-A y/«=0.25 (b) Section B-B j//«=0.50 (c) Section C-C y/»=0.75 (d) Section D-D y/«=0.95 

Figure 11. — Load distributions calculated for four streamwise sections of an untapered wing; m=0.6; 5/1 = 1.92. 
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The results of the calculations are presented in the form of 
values of P(Aplga) (equation (5)). 

The various components of lift are presented separately as 
calculated. In figures 11 (a) and (b),. the discontinuities in 
slope show the effect of the cancellation of the finite velocity 
Uq at the trailing-edge apex. The integrated part of the 
trailing-edge correction (component 2) has zero slope at the 
Mach line. The two outboard sections (figs. 11 (c) and (d)) 
are intersected by the Mach cone from the tip, as indicated 
by the finite drop in the load curves. Cancellation of the 
finite tip effect at the trailing edge (component 4 in both 
figures) results in a sharp discontinuity in pressure gradient 
along the reflected Mach line, at 91-percent chord at 
y/s=0.75 and at 78-percent chord when yls=0.9o. The 
cancellation of the trailing-edge corrections at the tip, which 
affects only the last section shown, results in another break 
in the load curve at 49-percent chord. Further corrections 
enter at the rear of the section as a result of successive can- 
cellations of the superposed pressures at the tip and trailing 
edge. Their effect has been only estimated. 

SWEPT-BACK WINGS WITH INTERACTING TRAILING 
AND LEADING EDGES 

When, as in figure 1 (b), the Mach cone from the trailing- 
edge apex includes a region ahead of the leading edge, the 
previously calculated trailing-edge corrections to u must be 
canceled in that region, since the}^ represent a discontinuity' 
in pressure which cannot be supported in the free stream. 
Thus there must be calculated a leading-edge correction, 
which is one of the previously defined secondary corrections. 
However, the location of the disturbed field ahead of the 
wing causes its influence on the wing to be so much more 
widespread than that of the other secondary corrections as 
to require more careful consideration. A new type of flow 
field is also required, as discussed in the following paragraphs. 

LEADING-EDGE CORRECTIONS 

Elementary solution for the region ahead of the leading 
edge. — In general, the elementary solution required for the 
cancellation of pressure in the plane of the wing ahead of the 
leading edge is one that : 

1. Provides constant streamwise velocity over an infinite 
sector bounded on one side by the leading edge of the wing 
(extended) and on the other by an arbitrary ray extending 
outward into the stream from some point yt, on the leading 
edge. (Seefig. 12.) 

2. Induces no vertical velocity, or downwash, on the wing. 

3. Induces no lift except on the wing and within the sector 
described in condition 1. 

At first glance these conditions would appear to be satisfied 
by the oblique solutions used at the trailing edge, if properly 
oriented with respect to the wing, and the same form of solu- 
tion might be expected to apply. In reference 4, however, it 
has been pointed out that the downwash connected with the 
latter solution remains constant over the wing only if the 
wing area does not include the line ?/=constant extending 
downstream from the apex of the element. In the case of 
the leading-edge element this condition is violated (fig. 12) 
and an additional term is needed to bring the downwash to 
zero throughout the area of the wing affected by the element. 



The solution applicable to this case has been given in refer- 
ence 4. The u component of the velocity in the plane of the 
wing is as follows : 


Ub r 

u=r. p. — cos 


_i (tg—m) (1 +tt) — {m — U) (1 +tg) 
(l+m) {t(,—tg) 


2m 


(l-f-m)f, 


(34) 


where Ui, is the constant streamwise perturbation velocity 
over the element, and U refers as before to a ray from its 
apex. The ray bounding the element originates at a point 
on the trailing edge and has been designated, from equation 
(24), as tg. When the correction is being made for the sym- 
metrical trailing-edge element, tg is replaced in equation (34) 

For brevity, the two parts of the correction function will 
be referred to as 


and 


C(U-r.p. cos ^ ’ 

The variation with U of these functions and the induced 
velocity (equation (34)) are illustrated in figure 12. 
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Leading-edge correction to the loading. — The single conical 
field of {Au)o will be considered first. (See fig. 13.) The 
velocity field to be superposed ahead of the leading edge to 
cancel the velocity (Au)o induced in the plane of the wing 
by the symmetrical solution (equation (20)) can be built up, 
as shown in figure 13, of overlapping constant-velocity sec- 



tors having one edge along the leading edge of the wing and 
one along the extended ray from the apex of the trailing 
edge. The magnitude of the constant velocit}^ on each ele- 
ment is dto or, from equation (20), 

(1*0 

/o^\ 

Applying equation (34) to the cancellation of the sym- 
metrical-correction velocities {Au)q ahead of the wing results 
in the following induced velocity increment at any point 
(a;, y) on the wing: 

(A 2 «)o=- [< 7 «o)+ 2 ?«o)] du ( 38 ) 

7T J TO dlo 

where tq is that value of for which — 1, and designates 
the most rearward leading-edge element containing the 
point Xy y within its Mach cone. In terms of x and t/. 


_ m{x+fiy) 
(a^+iS2/)-(l + m)Co 


(39) 


Integration of — dto is not feasible by elementary 

means. For graphical integration it is advisable to rewrite 
asrf(Ati)o to avoid the infinite value of the derivative 

at ^ 0 = 1 - 

The second term of the product in equation (38) can be 
integrated in closed form as follows: 


r ... _ -Am^^^UoK{k) I x + 0y 

Jtq dto ^ ° mt(l+m)K\mt)\ mx—fiy ^ ^ 

/ 2m,(l-To ) 


where 


and 


with 


,)(ro+m,) 




(41) 




0+rrit 


2 to 


The function Z(^, k) is tabulated in reference 16; a plot of 
Z{\(/ k) . 

k~syn~^ against ^ is given in figure 14. 

^Similarly, for each oblique trailing-edge element a (see 
fig. 15), a canceling field can be built up ahead of the leading 
edge by the superposition of sectors bounded by the leading 
edge and by rays ta from the apex x^, ya of the element a, 
and having a constant velocity of the magnitude 


^(^^)« 1 ^ (1 —a){ta-mt)-{mt—a){\ -Q 

bta (l-m,)(4-a) 


dt(i 
(42) 

(from equation (24)). If the symbol Aul.e. is used to 
designate the total leading-edge correction to the u com- 
ponent of velocity at any point, then the part due to canceling 
the field of a single oblique trailing-edge element a is 


dAu 


L.E. 


da 


da 


where 


=-f 


A(A«)„ 

d/. 


ma) + R{Q]dh 


m(?n,— g)(z+<3y)— m,Co(l +m)g 
(m, — g)(a: + (3i/) — m ,Co( 1 + m) 


(43) 


(44) 


is the value of for which t|,= — 1 and the leading-edge 
correction function vanishes. 

When the expression (equation (42)) for is sub- 

Old 

stituted in equation (43), it is again impractical to attempt 
to write a closed expression for the integral of the first term 
b(Au)a 
dta 

term is 


C{ta) of the product. The integral of the second 


x; 


d(Ag)a 

i>ta 

— 4m 


R{Qdta= 


Ua l{m, 

aV~ 


— g)(l— m)(a;-i-^y)— miCo(l-fm)(l— g) 


7t( 1 -)-m) 

K(ka) j^Vl+a 


Z{ia, kg) 

ka sin \ka 


(1 — m,)(mx— j3y) 

0 


^ ka sin \ko. 


X 


(45) 


976745—52 3 




14 


REPORT 1050 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


sin ** k 



40 50 ^ 

Argument , y, d€0 

Z(^,k) 

Figure 14.— The elliptic function — — • 
il: sin ^ 
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where 




m,){l — r„) 


-m,)(l + T„) 

- sin- >■» 

' \{ra-c 


■a){l+m,) 


'I'o 


■ -I /wi, 


(1 +Tfl) 


(1 +^«) 



Figure 15— Cancellation of the pressure field introduced ahead of the leading edge by a 
single oblique trailing-edge element. 


Then the total leading-edge correction to the velocity u at 
any point x, y is 


E, = (A2«)o + da (46) 

vhere (A 2 ti)o is obtained from equations (38) and (40), 
^ equations (43) and (45), and 


/_ (l+m)co— (1— m)(x + /3y) 

® (1 +m)mtCo— (1 —m) {x-\-0y)^^ 


(47) 


is the value of a at which ra{x, y, a) (equation (44)) is equal 
to 1. 

The last term in equation (46) will seldom be found to 
contribute any significant amount to the loading, but will be 
needed in calculating the leading-edge thrust. 


FURTHER CORRECTIONS 

Omitting for the moment any specification of tip location, 
it is in any case necessary, as seen in figure 12, to consider 
the effect of a further cancellation necessitated by the excess 


lift introduced behind the trailing edge by the leading-edge 
cancellation field. To compute this effect by the conical- 
flow method would be feasible only with the aid of high-speed 
computing machinery. The previously mentioned cancella- 
tion method of reference 5, being more direct, would be some- 
what easier to use in this connection, but the calculations 
would still be very lengthy. It will be shown by numerical 
example that the effect of the first cancellation at the trailing 
edge of the leading-edge correction, which is initially quite 
small, may be estimated with adequate accuracy when the 
section loading is considered as a whole, provided the fraction 
of the chord affected is not too large. 

If the product cot A is low or the aspect ratio high, still 
further cancellations will be required (see fig. 16) at both 
leading and trailing edges. It is clear that calculation of the 
effect of these further cancellations hy the conical-flow 
method is all but impossible. The doublet-distril)ution 
method of reference 5 does not appear to offer any consider- 
able advantage in this application since, in canceling lift 
ahead of a subsonic leading edge, it is necessary to find not 
only the pressure distribution to be canceled, but the asso- 
ciated sidewash distribution as well. 

It is apparent that an alternative method must be sought 
for describing the flow in the outboard regions of a high- 
aspect-ratio wing or a wing the sweep of which is large com- 
pared to the sweep of the Mach line. If the wing could be 
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extended indefinitely, it is known that the fiow must even- 
tually approach the two-dimensional subsonic flow, in accord- 
ance with simple sweep theory. The question then arises, 
can the flow at a distance of the order of a semispan from 
the apex of the swept-back wing be related to the two- 
dimensional asymptotic flow? While the flow field appears 
to be too complex to obtain an answer to this question on 
analytical grounds, numerical values, presented in the 
following paragraph, suggest a practical approach. 

NUMERICAL RESULTS (WITHOUT TIP EFFECT) 

Load distributions ha ve been calculated by the conical- 
flow method for three combinations of taper, sweep, and 
Mach number as follows: 

Untapered Tapered 

m= 0.2 0.4 0.4 

mt= 0.2 0.4 0.6 

These values of m and rrit represent, by virtue of the Prandtl- 
Glauert transformation, a variety of sweep angles at Mach 
numbers between 1 and 2; as for example, 0.2 would be the 
value of m for a wing with 63° sweep of the leading edge at 
a Mach number of 1.075, or 75° sweep at a Mach, number of 



(a) Section A-A /Sy/co=0.667 

Fioure 17.— !Load distributions calculated by the conical-flows method 


1.25. Similarly, m=0A would correspond to 45° of sweep 
at M=1.08, 60° at M=1.22, or 75° at M=1.80. The 
trailing-edge sweep angles at these latter Mach numbers, 
if m,=0.6, are 34°13', 49°, and 68°, respectively. 

Figure 17 presents the lift distributions at two stations 
of the tapered wing. Each component is plotted independ- 
ently in order to show" the magnitudes at the leading edge. 
Section A-A contains the intersection of the trailing-edge 
Mach line w’ith the leading edge, so that the value of the 
leading-edge correction is zero at the leading edge of this 
section. At points farther back along the leading edge, as 
at /3y/co=0.8, the correction is minus infinity. However, it 
is seen to increase to a small positive value within a fraction 
of the chord length at this station. 

At both stations it is necessary to estimate the effect of 
cancellation of the leading-edge correction at the trailing 
edge to satisfy the Kutta condition. Cancellation would be 
carried out by means of oblique elements of the type used 
previously (equation (24)) in canceling lift at the trailing 
edge. The pressure to be canceled is initially (i. e., at 
(%• 16)) zero. Then the lift induced on the wing 
by this cancellation ma}' be presumed to have the same 
general shape as the oblique trailing-edge correction of figure 



(b) Section B-B ^i//co=0.800 

for two streamwlse sections of a tapered swept-back wing; m=0.4; 7nt=0.6. 
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(a) Section A-A ^i//co=0.6T)7 (b) Section B-B /Jy/co=0.800 

Ficure 18.— Load distributions calculated by the conical-flows method for two stroamwise sections of an untapered winjr: m=0.4. 


11, falling along a modified inverse cosine curve from the 
value of the error at the trailing edge to zero, with zero 
slope, at the boundary of the region affected. With this 
boundary (the Mach line from the point X 2 , 2 / 2 ), it is possible 
to draw a satisfactory estimate (dotted curve) of the correc- 
tion needed to bring the pressure once more to zero at the 
trailing edge. 

The untapered wing with the same sweep (m=0.4) relative 
to the Mach lines is shown in figure 18, with the load dis- 
tributions calculated at the same stations. 

Four section lift distributions are presented (fig. 19) for 
3v 

m=0.2. At —=0.15 only the rear 60 percent is influenced 

by the subsonic trailing edge. The reflection of this influence 
at the leading edge alters the pressure over the rear 40 per- 
cent of the section. At section B-B, the leading- and 
trailing-edge interaction affects the entire section. A further 
reflection of this effect at the trailing edge must be estimated. 

At section C-C the influence of cancellation of the leading- 
edge correction at the trailing edge extends over the whole 
of the chord and any estimate of its magnitude would be 
necessarily arbitrary. Also, a second pair of reflections 
must be taken into account. The final pressure distribution 
has therefore been drawn as a band within which the true 


curve may be shown to lie. Its height is the error introduced 
at the trailing edge by the first leading-edge correction, 
except very near the leading edge, where an infinite negative 
correction is known to be introduced by the second leading- 
edge correction. The calculations were also carried out for 
i^y/co=0-45. The margin of uncertainty was found not to 
have increased by any appreciable amount. (See fig. 19 (d).) 

APPLICATION OF TWO-DIMENSIONAL FORMULAS TO CALCULATION OF 
LOAD DISTRIBUTION 

Correlation of two-dimensional and swept-back-wing 
loadings. — It is apparent from the calculated results that, 
whenever the plan form and the Mach number are such that 
the trailing-edge Mach line intersects the leading edge, the 
load distribution behind the Mach lines from the point of 
intersection resembles in shape the theoretical load distribu- 
tion over an infinitely long flat plate in incompressible flow. 
However, as the results have been plotted, the quantitative 
agreement is not good, particularly in the case of the tapered 
wing. On the other hand, if the load distributions in 
cross sections normal to the stream are examined, a near 
proportionality of the curves is observed. In order to 
determine the factor of proportionality, it is only necessary 
to find the ratio of the strengths of the singularities at the 
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(a) Section A-A /9.v/':o=0.15 

(b) Section B-B fiylco-f).25 

(c) Section C-C ^.v/fo=0.35 

(d) Section D-D 0ijlco=O.ib 

Fioure 19.— Load distributions calculated by th^ conical-flows method for four streamwise sections of an untapered wing: 7n=0.2. 
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leading edge. T\v^n an approximate expression for tlie 
loading on the outer portions of a high-aspect-ratio wing 
can be obtained by adjusting the two-dimensional loading 
by that factor. 

Both the swept-back-wing and the subsonic two-dimen- 
sional loadings approach infinity as the reciprocal of the 
square' root of the distance to the leading edge. In sections 
normal to the stream, the distance from an\^ point x,y to the 

leading (*dg(‘ may be written^ The value* at the 

leading edge of the coefficient of {mx—fiy)~^^‘^ will be referreel 
to as the strength of the leading-edge singularity. 

The subsonic two-dimensional perturbation velocity has 
the form 


“-"W 


(48) 


where rj is the distance to the leading edge, expressed as a 
fraction of the chord, and B is a constant. If the section of 
the swept-back wing is taken perpendicular to the stream 
{x constant), the chord length is 


and 


- [mx — nitix — Co)] 




mx — mt{x^Co) 

Substitution for rj in equation (48) gives 




0y—m,{x—co) 


mx — Py 

Then the strength of the leading-edge singularity in u is 


(49) 


(50) 


(51) 


B y'm x — m,(x—Co) 


(52) 


and, for each value of a from 0 to that value ao' which makes 
Ta equal to one, 


(^^C , -4m 

aa = 

(ta 


u„ /(m, — a)z — (1 — a)m,Co , ,, 

a V X 

r j n . . 

L"' ka sin ypa ^ ka sin 


wh(*r(‘ To (equation (39)) and Ta (equation (44)) reduce to 

mx 

To = 

and 


X — Cq 

{mt — a )mx--mtCQa 
{mt — a)x — mtCQ 


and the arguments and moduli of the elliptic integrals 
follow as for equations (40) and (45). 

The coeflficient of {mx — at the leading edge is, 
therefore, in region II, figure 16, 

+ (56) 

with «o' reducing to 

// m \ CQ—{\—m)x 

a. (»,yx)=m, (57) 

Equating the two coefficients, expressions (56) and (52), 
gives for any one section 

B= — ^7 r r^' ( 58 ) 

^mx — m,(x—Co)\_ Jo da J ^ ' 

For convenience, a nondimensional cofficient 


The leading-edge singularity in the loading on the swept- 
back wing is initially (region I, fig. 16) that in the triangular- 
wing loading. Introduction of the leading-edge corrections 
to the load, in region II, reduces the strength of the singular- 
ity there through the terms /^(^o) and i?(0- (The inverse- 
cosine function is always finite.) The coefficient of 
{m 2 — in Ua is, from equation (6), 


mxup 
^Jmx + Py 

reducing to 

C^=Uoyl~ (53) 

at the leading edge. 

From equations (40) and (45), decrements to this coeffi- 
cient maybe derived for the portion of the leading edge just 
behind the intersection Xi,yi with the trailing-edge Mach 
line, as follows: 


(AOo^ 


— ^mupK{k) I 

l^^'(m,) V 


mx 

1 +m 


Z{4^M 


(54) 


is defined, so that 


B=Va 


(t{x)^j'co 

^Jmx — mt{x—c^ 


By substituting for B in equation (51), the loading on the 
outer portions of a swept-back wing is obtained as 


u ^ / Co[fiy-mt{x-Co)] 

Va, \ [m.x — mt{x—Co)]{mx — Py) 


(60) 


Numerical results. — The closeness with which the fore- 
going procedure predicts the theoretical loading over swept- 
back wings is indicated by figures 20, 21, and 22, where the 
previously calculated load distributions are compared with 
those calculated by equation (60). Even in the case of the 
highly tapered wing, the agreement is seen to be good. 
At the most inboard section of the m = 0.2 wing (fig. 19 (a)) 
there is, of course, no agreement over that portion, forward 
of the 60-percent-chord point, where the flow is essentially 
conical. At station B-B, however, the agreement is very 
good. At sections C-C and D-D, where the exact theoretical 
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(a) Section A-A 
Cb) Section B-B 

Figure 20— Load distributions on the tapered wing as calculated by the 
conical-flows method, compa*ed with the two-dimensional approximation. 


loading had not been determined, the two-dimensional-tvpo 
loading lies within the band prescribed b}' the conical-flow 
calculations. Since the discrepancy between the corrected 
two-dimensional loading and the exact theoretical distribu- 
tion is already, at section B-B (fig. 22 (b)), less than the width 
of the bands in figures 22 (c) and (d) and must diminish to 
zero at infinity, it may be supposed that the corrected two- 
dimensional curve is at least as satisfactory an approximation 
to the correct curves at sections C-C and outboard as 
at section B-B. It is probably more satisfactory than can 
be obtained by a limited application of the conical-flow 
method. 

The load distributions derived by simple sweep theory 
are included in the last part of each figure to show the magni- 
tude of the plan-form effect and also, in the case of the un- 
tapered wings, the curves that the load distributions must 




(a) SecLiJii A-A 

(b) Section B-B 

Figure 21.— Load distributions on the untapered wing, 7n=0.4, as calculated 
by the conical-flows method, compared with the two-dimensional approxi- 
mation. 
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(a) Section A-A 

(b) Section B-B 

(c) Section C-C 

(d) Section D-D 

Figure 22.— Load distributions on the untapered wing, 7n=0.2, as calculated by the conical-flows method, compared with the two-dimensional approximation. 


approach as the distance from the plane of symmetry is 
increased. In figures 21 (b) and 22 (b), comparison is also 
made with results of the slender-wing theory of reference 3. 

Discussion of the a function. — In the calculation of the 
pressure coefficient at points toward the rear of most of the 
sections considered in figures 20, 21, and 22, it was necessary 
to find (t{x) for values of x greater than Xz (fig. 16). In 
deriving a-(x), it was mentioned that expression (56) applied 
to region II. In region III, the strength of the leading-edge 


singularity is affected by further modifications of the flow 
taking place in region Ilb, so that additional terms in a(x) 
should be considered when x is greater than X3. Evaluation 
of these terms by presently known methods would require, 
as suggested earlier, the aid of high-speed computing ma- 
chinery. However, the successive terms are all initially 
zero and enter with zero slope at X3, zero slope and curvature 
at X5, and so on, so that the three-term expression for <r given 
by equation (59) may be used with satisfactory accuracy 
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for some distance beyond the last value of x for which it is 
strictly valid. In practice, the third term in equation (58) 
may also be neglected for values of x only slightly greater 
than Xx. 


Charts have been prepared (fig. 23) giving 



O' 


as a function of ^ — — for several values of the ratio m/m^. 
^0 

This last parameter is the ratio of the tangents of the semi- 
apex angles of the leading and trailing edges and is constant 


through the Mach number range for any one wing. The 
value of Xi is readily determined: 




Cq 

l — m 


(61) 


The curves were computed using equation (59) and are 
therefore exact only up to x = xz (shown by a vertical mark 
on each curve). Cross marks are drawn at the points x=X 5 
to indicate a more practical limit to which use of the curves 



O .4- .8 /.2 f.6, 2,0 2,4 2,8 3.2 oo 

sc-x^ 


Co 

(a) mlmt=l.O 



(b) mlmt=0.9 

Ficure 23.— Charts for determining <r, the strength of the leading-edge singularity. 
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(cj m/mt=0.8 


(d) 77»/mi=0.7 



(e) m/mt=0.(\ 


([) mlmt =0.5 
Figure 23— Continued. 


(g) m/mt**0A 


may be extended. (These points are off the scale for mt=0.S 
and 0.9 in figure 23 (a).) When the wings are untapered 
(m/m, = 1.0), asymptotes 


Il — m 

1 

/ m 

Vi 


derived from simple sweep theory, may be drawn. 

The curves, for the most part, are regular enough to permit 
interpolation within intervals of 0.2 in m,. However, at 
/n,= 1.0 the lines diminish to a point on the vertical axis; a 
curve for m,=0.9 was therefore inserted in the charts for 
values of m/m, equal to or greater than 0.5. When m/m, is 
less than 0.5, m = 0.9 represents, if the leading edge extends 
beyond x^yi^ such extreme taper that the successive reflection 
of the Mach lines (at xs, Xs, . . .) take place within a very 
small fraction of a chord length and no useful curve can be 
drawn. No curves are drawn for values of m, smaller than 
0.2 because of the tip-interference limitation mentioned in 
the introduction. 

Calculation of tip effect. — The foregoing assumption of 
two-dimensional flow can be extended to give fairly simple 


1 approximate formulas for the tip effect on a high-aspect- 
ratio wing. It is assumed that the velocity distribution to 
be canceled in the stream outboard of the tip is cylindrical; 
that is, is an extension of the velocity distribution calculated 
for the tip section along lines parallel to the leading edge. 
For this purpose the approximate load distribution given 
by equation (60) is used, still further simplified by assuming 
a to remain constant at its value at the leading edge of the 
tip section. (Where the wing is tapering to a point and a 
is changing very rapidly, the tip region is so small that the 
entire calculation of tip effects could probably be omitted.) 

The assumption of constant a results in a failure to cancel 
exactly the lift along the tip. The assumption of cylindrical 
flow, while reasonable for the untapered wing (compare fig. 
21 (a) with fig. 21 (b), for example) would appear to be too 
drastic for the tapered wing, where neither the chord nor 
the loading remains constant. However, as has been men- 
tioned earlier, the major part of the tip effect results from 
the cancellation of the infinite pressure along the leading 
edge, and this part will be accurately calculated. The 
effect of the residual lift on the rearward portion of the tip 
section and in the stream should be small. 
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The distribution of perturbation velocity at the tip station 
y=Sj with the simplification of constant cr, is, from equation 
(60), approximately 


• , V T 7 / [fis — rrit (Xc—Cq)] Cq 

U (Xc,s) (Ts (^Xc—Co)] (mXc—fis) 


(62) 


where Xc^s are the coordinates of a point on the tip and a, 
is the value of <r at the heading edge of the tip section. 

This expression may be more conveniently written in 
terms of the parameter 


= 1 -— 
^ nit 


and the variable 


Xc— 


^c=- 


m 


(63) 


(64) 


which is the distance of Xc,s from the leading edge (see fig. 
24) expressed as a fraction of the tip chord Ct. Since 

V rrit/ m 

equation (62) may be written 


u(?c,s)== 


.Fa 


^[m\ 


I 1-gc 

V(l-M^.)^c 


(65) 


where X is the taper ratio Ct/Co- 

If the velocity distribution u is assumed to be constant 
beyond y=s along lines parallel to the leading edge, it may 
be canceled by the superposition of conical flow fields of 
which the constanr-velocity regions have one edge along the 
tip and the other parallel to the leading edge, with apexes 
displaced along the tip Dy increments in The velocity 
induced at a point x,y by each such eleme 
(equation (11)) 


where 


and Uc is the velocity on each sector. 

Following the procedure used in deriving equation (14), 
the corresponding equation may be written for the pres- 
sures induced by canceling the cylindrical flow 


(^) = 


tip 


— Vm (1 + ffl) (a: — X q) 
irVa 


r- 

Jfis I 


U{XcyS) dXc 


[(l + m)x—Xo—mXc\ -ylxQ—Xc 


( 68 ) 


where Xq,s is the intersection of the Mach forecone from x,y 
with the tip. 



If the distances of x,?/ and Xq^s back of the leading edge, 
measured as fractions of the tip chord, are 


x.y by each such element 

would be 



(69) 

Ue m^-tc-\-2mtc 

— cos ‘ 

IT tc—m 

(66) 

and 


(70) 

f fiiy-s) 

"" X — Xc 

(67) 

it can be shown that 

(1+m) (x— Xo) = mc,(^— Jo) 

(71) 


from which equation (68) can be written (with the substi 
tution for u{Xc,s) from equation (65)) 


(rj 

1 _ 1 

r*> di, 

1 1-ic 

tip ^^^|n^\ J 

1 

o 

(fo-{.)(l-M{c)«« 


In integrating equation (72), three cases must be distin- 
guished: (1) (always true for the untapered wing), 

(2) (when the point lies more than a tip- 

chord length behind the leading edge), and (3) “ 
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(a possible condition for some points near the trailing edge 
of a highly swept or tap(*red wing). 

In the first eas(‘ (^<Cl) 


a/np 1 — " 


where 


-= 


^o) 




(73a) 


Mf) 


and Ao is the function (equation (16)) plotted in figure 6. 


In the second case 




equation (73). On the trailing edge of an untapered wing, 
^=1 and 



There is no corresponding simplification for the tapered wing. 

Numerical examples, tip effect.— Equations (73) and (31) 
have been used to calculate the tip effect in two cases, namely: 
m = m,=0.4, Ps=0.94Co; and 7^ = 0. 4, m^=0.6, ^s=0MCo. 
ddie tip effect has been calculated for each wing at Py=0.8co, 
where the loading was previously calculated (figs. (17b) and 
(18b)) assuming the wing to extend indefinitely. The tip 
locations were selected so that in each case only one reflection 
of the primary tip effect affected the section at 0y=O.8Co. 

Figure 25 shows the results of the calculations. The heavy 
solid curve in each case was calculated entirely by the 
corrected two-dimensional theor}' — that is, by equations 
(60), (73), and (31). As a check on the accuracy of the 
cylindrical-flow approximation for the flow outboard of the 
tip location, the accurate theoretical loading was calculated 


f ^0 P /7 t” 7t iTT ^) ~] 


(73b) 


where 


/ • -1 l^t 

^ 2 = sin — 


(1-M« 




and Z is the function (equation (41)) plotted in figure 14. 
In the third case 

\'^ M/ 

— (^s Ao(^ 3 , IC) 




V»iX? sin \{^3 


(73c) 


where 




Along the Mach line from the leading-edge tip all three 
equations reduce to the value 


Au* — cr. 


(74) 


By the procedure just described, approximate cancellation 
of all pressure differences outboard of the tip has been 
effected, but the pressures induced by- such cancellation 
now violate the condition of zero lift in the wake. Approxi- 
mate cancellation of the induced pressure differences in the 
wake region can be accomplished, as before, by making use 
of the known value of the tip-induced velocity at the trailing 
edge of the wing, but assuming the entire error to originate 
at the leading edge of the tip. Equation (31) is directly 
applicable, with At/* given by equation (74) and {Au)tiphy 



(a) Untapered wing; m=0.4; fiA = l.BS. Section at fiy=0.8co, or 85 percent semispan 


Fioure 25. — Load distribution over streamwisc section near lip as calculated by two-dimen- 
sional formulas, compared with more accurate theoretical values. 
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(b) Tapered wing; m=0A; m <=0.6. Section at fiu=0.8co, or 93 percent semispan 
Figu re 25— Continued 


for one point within the region of influence of the tip in each 
case. The procedure employed for the exact calculation was 
as follows: 

The accurate loadings with no side-edge effects had already 
been calculated, as has been noted, by the conical-flow 
method. A primary tip correction was calculated for each 
case by equation (15). This correction is the effect of cancel- 
ing the unmodified triangular-wing loading off the tip sta- 
tion. The remaining pressure differences to be canceled 
consisted of those introduced by the leading-edge and trail- 
ing-edge corrections. These pressures were computed by 
means of equations (2(3) and (46) of the present report and 
canceled by the method of reference 5. 


The results are designated by the circled points on each 
figure. At the point at which the section enters the tip Mach 
cone in each case, a second circled point indicates the accurate 
theoretical loading. The value differs from that calculated 
by the approximate formulas only as the two loadings without 
tip effects differ. 

It may be pointed out in concluding this section on load 
calculations that, while the formulas have been developed for 
plan forms with streamwise tips, the procedure may be 
adapted by obvious means to raked tips as well. However, 
in every case the deviation in the tip regions of the physical 
flow from the assumed potential flow must be borne in mind. 

Ill— LIFT ® 

GENERAL PROCEDURE FOR CONICAL FLOWS 

The total lift for any wing is, of course, the integral of the 
loading over the wing area. In general, however, it is difficult 
to obtain an analytic expression for the lift by a direct 
integration of the lift distribution. In the conical-flow 
method, advantage may be taken of the simplicity of the 
component fields by integrating the lift associated with each 
one and then combining the results in the same way as the 
pressure fields. 

Conical elements of area are employed for the integrations. 
These are infinitesimal triangles bounded by two adjacent 
rays of the conical field and the intercepted boundary of the 
wing plan form. Over each of the infinitesimal triangles the 
velocity u of the conical field will be constant. Thus it 
remains only to perform a single integration, with respect to 
the conical variable of the field, to obtain the total lift 
associated with that field. 

GENERAL FORMULA FOR THE LIFT INDUCED BY A SINGLE TIP ELEMENT 

The lift (AL)a induced on the wing by a single canceling 
tip element is obtained first. Although the notation of the 
solution (equation (11)) used to cancel the triangular-wing 
loading is employed, the derivation will hold generally for any 
canceling element bounded on one side by the tip of a swept- 
back wing, since no use is made of the fact that the other 
boundary of the element passes through the origin of the 
x,y axes. We write 

(AL)«= 2 P F J_"^ (A«)„ ^ dta (7 6) 

where (Aw)„ (equation (11)) is the streamwise increment of 

velocity induced by the canceling field and ^ dta (fig- 26) is 

the element of wing area S for integration. For simplicity it 
will be specified that the Mach cone from the apex of the 
element does not include the apex of the trailing edge nor any 
part of the opposite tip. Then (see fig. 26) 

dS / x, — Xg \^ (77) 

dta 2 ^ \mt—ta) 


• « It may be noted that, as a result of the reversibility property (reference 17), the formulas 
for the lift given herein for swept-back wings are equally applicable to the swept-forward 
wings having the same plan forms but reversed in heading. 
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Fioure 26.— Sketch for the determination of lift induced by a single tip element. 


Substituting from equations (1 1) and (77) and integrating 
by parts, we obtain 

Ua{xt — Xafg{a) (78) 

where 


(j{a)= 




a(l ~\~d} 



(79) 


and x,—Xa is the distance of the apex of the element from the 
t railing-edge tip. 


GENERAL FORMULA FOR LIFT INDUCED BY OBLIQUE TRAILING- EDGE 

ELEMENT 

With the notation of equation (24) for the velocity field of 
an oblique trailing-edge element, and on the assumption that 
the Mach lines from the apex of the element do not cross the 
leading edge, the formula for the elementary area of integra- 
tion with apex on the trailing edge (fig. 27) is written 


2t 


did 


(80) 


where s— is the spanwise distance from the apex of the 
element to the wing tip. Then the lift associated with the 
element is 


rs nT/W-r (l-a)(f„-m,)-(m<-a)(l-Q^ , 

(A/A=2pF-J^_cos ' dt/^- 


Integration of equation (81) gives 


(81) 


{AL)d=p Vis-ydYfi 


r /(W(— g)(l— g) 

g LV m, 




WING WITH SUBSONIC LEADING EDGE 

UNCORRECTED LIFT 

First, the uncorrected triangular-wing loading (equation 
(6)) is integrated over the wing plan form. The element of 
area is a triangle formed by two rays from the leading-edge 
apex g and a-\-da and either (1) the trailing edge of the wing 
or (2) the wing tip, accordingly as a is less than or greater 
than at, the value of a corresponding to the ray through the 
tip of the trading edge. (See fig. 2.) In the first case the 
differential of area is 


mtW 


and in the second 


2^{vi, — ay 


da 


2a‘ 


da 


so that the total uncorrected lift is 

r_4pFf C°‘ . 1 , 

From the geometry of the wing, the relation 

m,Co=^ (m,-a,) (84) 

dt 

may be deduced. With this substitution, and the substitu- 
tion for Ua from equation (6), equation (83) may be inte- 
grated to obtain 


Lo_ 4s* /3«o \ 

[m*(m,— a,)*| 

r 1 

( -1 ny-rritat 

q;a may Va j 

1 my—m^ 1 

\_-\/my—m^ 

7n(m,-g,) 


cos’ 


w,/ m,j mt—my ) 


(85a) 


28 


REPORT 1050 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


When mt=mj this reduces to 
qa Smat Va 

It should be noted that, for a given plan form, 
with Mach number only as Uq. 


(85b) 

varies 


WING WITH SUPERSONK/ TRAILING EDGE (TIP CORRECTION) 

Proceeding to the calculation of the tip correction to the 
lift, we integrate the change in lift (AL)a (equation (78)) 
induced by each element a over the range ai<a<m. The 
quantity u^(m) is substituted for Ua of the initial canceling 

element and da for Ua for the remaining ones. As in 

calculating the tip-induccsd pressure correction, the difficulty 
is encountered that u^{m) is infinite, and therefore the total 
lift correction must be written in terms of limiting values. 
Following the substitution 

X, — Xa=^s(^ — i) (86) 

\at a/ 

in equation (78), it is convenient to define the function 


G{a) 


(a — a,y 


g(a) 


Then the total induced lift may be written 
AL=2pV'm^* /3s* lim 1^— «4(a)(?(a)4- 

Integrating by parts results in cancellation of the first 
term inside the brackets Since (7(a,) is zero, equation (87) 
reduces to 

(—^ = - 4m,* /3s* r G\a)da (88a) 

\qot/tip J^t 


where 


G'(a) 


a—Qt 

a^a^ 



a — at 


2(/n , - o) A (m, + m,*)(a + a*), 


] 


is the derivative of G{a). 

Equation (88a) has been integrated (appendix C) in 
terms of an incomplete elliptic integral of the third kind. 
If the necessary tables are not available, it may be prefer- 
able to integrate numerically.^ In that case it is noted that 


1 (I .y a 

, sm ^ — 

c/a m 


and da is rewritten as 



Equation (88a) then becomes 

Ih this way infinite values in the integrand are avoided. 

WING WITH SUBSONIC TRAILING EDGE 

The expressions (equations (85)) for the uncorrected lift 
apply regardless of whether the trailing edge is subsonic or 
supersonic. The formulas for the tip correction may serve 
as a first approximation when the trailing edge is subsonic 
if the accuracy of a second correction is not required. For 
that purpose the special value for the untapered wing will 
be of interest: 

If the wing is untapered the elliptic integrals in equation 
(88a) (see appendix C) reduce to the first and second kind 
and the primary induced lift may be written in the following 
closed form: 


/Ai\ 4^^ PUq 

\qccjtip~~ ^mat^ Va 

2m^{l-\-aty m-\-at 

a/(l+m)^(m + ad 2 

where 

. I / 'iTi at I 1 / 1 'oi 

\L=sm~\—7 ^. — X and k = ^ — ^ — 

^ \m(l+a«) V 2 

The primary tip correc tion, however, is usually quite large. 
It may therefore be desirable to take into account the sec- 
ondary correction resulting from its cancellation at the sub- 
sonic trailing edge. Rather than compute a single secondary 
correction to the lift, as an additional item, it is again found 
advantageous to treat each superposed field individually, that 
is, to cancel each conical tip field at the trailing edge and 
find the net effect on the lift, then integrate over all the 
tip elements for a combined primary and secondary tip 
correction. 


{m — aj)^ r (1 — ~| 

m(l +mY \_at 2(m + a,)J 

2{m + at) at (m — a,)^ “| ir/, 

a, m 2ma,(l+m)J My 

(88c) 

Tip correction with subsonic trailing edge. — For the can- 
cellation at the trailing edge of a pressure field originating 
at a point Xa, s on the wing tip, equation (82) is applied, 
with the parameter a, which defines one boundary of the 
oblique canceling field, replaced by referring to a ray from 
XajS. The velocity Ua is the gradient 


of the field (equation (11)) to be canceled. The distance 
from the apex of the canceling field to the wing tip is express- 
ible as 

Xt —Xg TYlttg 

P rflt — tg 




’ Or see reference 10 for an approxiioate formula valid when m is close to one. 
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Then tlie effect on the lift of canceling the single field from 










P U 

dta 


(89) 


m, w, J " 

which may be integrated, after substitution for to give 

dlQ 

&irnt—a) nii 


^ (1 +a) 




in which 
and 


(1 +mt) 

— (90) 
7rit L sin ^ J) 


\l+m, 

, ■ -1 I m, — a 

f=sm \ — — 

V w,(l — < 


If the foregoing result is added to the lift associated 
with the original tip element, given by equations (78) and 
(79), it is found that the latter lift is exactly canceled by the 
algebraic terms in the reflected lift, leaving 

(AA=p Vm.^Ua r 

/3(m,— a) |_ V w,(l + m,) m, sm ^ J 

(91) 

for the lift induced by one tip element and its cancellation 
at the trailing edge. It will generally be found that further 
steps in the cancellation process are unnecessary- for engineer- 
ing accuracy. 

For the total tip-induced correction to the lift, it is neces- 
sary to write as before 

{^L)nv=2pVm,'‘0sHim^—UJ^{a)J{a)+j ^J(a)daj (92) 
where J{a) is 

(a-a,y 1 |~ / a(l+g) ^ a Ap(,f^,Ar)1 

a,V + m, sin if- J 

An integration by parts reduces equation (92) to 

f*m 

(AL)up=-2p u^{a)J\a)da (93) 

with 

di d mt—d i\_\d rrit — d/ 

{ 1 _ l + m, \~| / a(l-hq) ^ ,, , 

2(ld-o)\g ' 

1 —d^ V rnt{l+mt) ' 

A f A/ I m,-aA (l-m,)(a-ad "] Ap(^, k)\ 

rnti\d~^ Mt — d 2(1— d)(mt--d)J sin^ ) ^ ' 


If the wing is un tapered, J'(d) becomes indeterminate when 
d = m. The limiting value is 

(S+i5S)- 

Further integration must be performed numerically. In 
order to avoid infinite values in the integrand, note again that 


/ \ d • \ d 

u^(d)=mu^ — sin~^ — 
dd m 


(95) 


so that equation (93) may be rewritten 

m 

Trailing-edge corrections. — In deriving the trailing-edge 
corrections to the total lift, primary and secondary effects 
will again be combined. Further corrections will be 
omitted. 

For the symmetrical wake correction, the element of area 
is obtained from equation (80) by setting ?/« equal to zero, 
and substituting % for ta. Then the decrement in lift in- 
duced by the application of the symmetrical canceling ele- 
ment at the trailing edge is, from equation (20), 


or 


Wo r 

mi 

(97) 


-4s' r 

1 ^ 1 

(98) 

Ml V a\_ 



where 


7T 


The effect of canceling the pressure field induced by the 
symmetrical wake correction at the wing tips is obtained 
with the aid of the previously derived formula (equation 
(78)) for the lift associated with a single tip element. The 
parameter defining the boundary of the canceling tip ele- 
ment is now ^0 instead of a, and the velocity on the canceling 
sector is 


d (Au)q 
dto 


dto — 


Uq dto 


K'(mt) ^J(l—U^ (to^—rrit^ 


(99) 


The distance from the apex of the canceling sector to the 
trailing edge may be expressed as 




(---) 

\mt tj 


( 100 ) 
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SO that the secondary effect of the symmetrical trailing-adge 
correction becomes 

{A^L)o=2pVm,'^fis^ r — rY /■ — — 

(101) 


or 


(A2L)o_4s 2 |3j^0| j 


qa 


with 


;[ 


m, Va\_ Ko (m,) 


2^/^ Ko(A:)-£;„(m ..... 

Ko'im,) J 


k-yj'- 


— rrit 


Addition of this seccndary correction to the primary 
effect given in equation (98) results in the single correction 

/AL\ ^ — 8^^ Puq I 2 Ko{k) — Eo{k) 

\qa/o Va^l+rrit Kq (rrit) 

By a similar procedure, the effect of canceling one of the 
oblique trailing-edge fields at the tip is readily obtained and 
added to the primary effect given by equation (82) to yield 

{M\=pVp (s-VaY ^ r Eo{k) V(/n,-a)(l-a)- 

a Vi'WjL 

rrit — a ll-^a "1 

with k = , ^ -and ^==sin~^ ^ ^ the combined 

^ l+rrit ^ \mt(l+a) 

primary and secondary (correction to the lift due to a single 

oblique trailing-edge cau.cellation. 

For the total correction to the lift due to cancellation of 

the gradient of the triangular-wing loading in the wake, 

equation (104) is integrated graphically or numerically 

across the span as follows: 

AL —4mpuo (s—yaY ri;- /7 \ 

2«“ VaJo 

rrit—a h-\-a "1 

(105) 

Numerical examples to be presented will show this com- 
ponent of the lift to be very small, in general. 

WING WITH INTERACTING LEADING AND TRAILING EDGES 

In computing the load distribution it was foimd that, 
when interaction takes place between the flow fields of the 
leading and trailing edges, the wing plan form appears to 
comprise two principal ]*egions separated (see fig. 16) by the 
Mach line arising at the point of intersection Xi^yi of the 
trailing-edge Mach line* and the leading edge. Ahead of 
this line (region I) th(c flow is most readily described in 
terms of conical fields. Behind this line the flow is more 
nearly two-dimensional. On this basis, the total lift will 
be found in two parts, using for region I the conical-flow 


expressions for the loading, and for the remainder of the 
wing the quasi-two-dimensional approximation. 

LIFT ON INBOARD PORTION OF WING 

The uncorrected triangular-wing loading will first be 
integrated over region I, shovm shaded in figure 28. For 



Figure 28 .— Inboard portion (region I) of high-aspect-ratio wing. 

this purpose the region is considered in two parts, separated 
by the ray ch. from the wing apex to the point 2 : 2 , 2 / 2 - When 
a is less than 02 , the element of area is as before 


mYcY 


When d>a 2 , the element of area is 


da 


(IH-m)V . 

2)3(1 -m)2(l+a)2 

Thus, the uncorrected lift in the entire shaded region is 

( 106 ) 
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or 


Lo 4m Uo 2 ( 1+^ r 1 
qa- p j(l-mfLVT3 

rrit^ / m 

m,^— m^\mt mi — ai ) 

{ -gj— yrcos- 

mv L 


cos”'- 


m^ + g2 Vm^— a2^ 
m(l+a2) 1+C2 


+ 


with 


d2- 


+ 

m^—mia2 
m(m, — tt 2 ) 

2mmt 


-cos 


■'3 


1 +m + m^ — m rrii 


(107) 


(108) 


When rrit—m (untapered wing), the second part of equation 
(107) becomes indeterminate. In this case, 


u 

qa 


4m Uo 2 ( 1 f 1 


^ m^+a 2 Vm^— a2^ 
1 . 


— cos . 

m(l+a2) 


+ 


(109) 


The trailiiig-edge corrections to the loading are to be 
integrated over the part of the shaded region behind tlie 
trailing-edge Mach lines. Integration of the symmetrical 
wake correction (equation (20)) yields 


(AX)o. 


— IGm^Co^ 


Uo 




qa |3(1 + mO(l — m)^ Va\ 

For each oblique element, the reduction in lift is given by 

1 dAL , 

j — aa = 

qa da 


/SFal+a^”"' LV 2(2^a) ^ 


with Ua=^^da, 


rritCo 

rrit—a 


and 


/1+m . \ 


^0 


+ rrit 


The total lift in region I is then given by 


(-)■ 

\2«/l 


Lq , (AZ)o_| 2 


-f 

a Jo 


^*2 dAL 
da 


da 


( 111 ) 


( 112 ) 


(113) 


qa qa qc 

The quantity is plotted against m^ in figure 29 for 

several values of the ratio ml nit, 

LIFT ON OUTER PORTIONS OF WING 

In order to find the total lift (except for tip losses) on the 
remainder of the wing (fig. 30), a double integration with 
respect to x and y is performed on equation (60). A first 
integration, with respect to i/, yields for the indefinite 
integral 

r ^ f7\ / (mx-^y)(m,Co-m<x+'j^ , 

J Fa P ^ IV mtCo-{mt-m)x + 

^mtCo—{mt — m)xian-\ — 

V m,Co— m,x + /32/J 


The values of fiy to be substituted as limits in equation 
(114) are indicated in figure 30. Along the leading edge, 
the right-hand member of equation (114) reduces to zero; 
along the trailing edge it becomes 

2 Vco ym^Co— (m, — m)x 

Then the total lift on the outboard region (both wing 
halves), except for tip losses, is 

f£' ,(x)j.dx-£ 

tn 

(115) 

where 

/j = ^'mtCo— (mt — m)x 

/2=V(1+«^)(2: — a^i) Jt=^/m{x—fislm) 

j3=-yl{i+m,){x2—x) /s= Vm,(x,— x) 

The indicated integrations may be performed numerically 
or graphically, using values of <t(x) taken from the charts of 
figure 23. 


TIP-INDUCED CORRECTION TO THE LIFT 

In deriving a tip correction to the lift, the same simplifying 
assumption of completely cylindrical flow will be adopted 
concerning the pressure field to be canceled as was used in 
obtaining a tip correction to the loading. As in the preceding 
section, a combined primary and secondary tip correction 
will be derived. All further corrections will be omitted. 

If the notation of equation (63) is used the distance from 
the apex Xc,5 of a canceling element to the trailing-edge tip is 
Ct(l — fc), and the lift induced by the element and its cancella- 
tion at the trailing edge is, from equation (91), 


(AL)c=pVmt^Uc 


fi(mt—m) 


V m{\ ~t~ m) rp / T \ ^ A.q(^, At) "] 

m,(l + mt) ® m, sin i/' J 


where 


as before, and 


Vl+m, 


, . / mt — m 

^=sm \ — -7Z — 

\m,(l — 


m) 


(116) 


since the outer boundary of each element now has the slope 
m 

T 

It is seen that only Uc and (1 — ^c)^ in the coeflScient of 
equation (116) vary with the element. For the first element 
(f^=0), the velocity Uc is the initial value of the uncorrected 
velocity along the tip section given in equation (65), and, for 
the other elements, the differential of that velocity. Then 
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Figure 29.— Chart for the computation of lift in region I. 
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Figure 30.— Boundaries of outboard region of high-aspeci-ralio wing, for use as limits of 
integration in equation (114), 


Figure 31.— Chart for the correction of the lift for tip effect, using two-dimensional 
formulas. 


the combined primary and secondary tip correction on both 
wing halves is 




r l m{l+m) 

|8 (;n, — m) -v'mX L V wi<(l + 


( 117 ) 


Integration by parts gives, finally, 

AL 8vm,^Ct^-</rn F /m(l+m) 

L V m,{l +m,) 

(118a) 


If the wing is untapered, equation (118a) takes on the value 


Ya=~W^) \‘^rtiK,(k) + {\-Zm)E,{k)\ (118b) 


Except for the occurrence of a„Ct and X in the coefficients, the 
tip correction obtained in the foregoing way is a function of m 
and rrit only, independent of the tip location. Values of 


have been plotted in figure 31 in a form similar 
to the chart of (%• 29). 


APPLICATION OF LIFT FORMULAS 

CASES COMPUTED 

The lift-curve slope has been calculated for two families 
of un tapered wings with varying aspect ratios as follows: 


m = 

=0,2 

m= 

= 0.4 

fi - 


fi - 


Co 


Co 


0.3 

0.6 

0.6 

1.2 

.4 

.8 

.8 

1.6 

.6 

1.2 

1.2 

2.4 



1.6 

3.2 


and for two tapered wings: 


m==0.4, OT/=0.6 

fi- 

0A 

X 

co 



0.6 

1.6 


.8 

2.4 



It should be noted that the untapered-wing cases (except for 
the last one under m—0.4) represent three wings of fixed 
geometry at two different Mach numbers such that is 
doublei^ in going from the first to the second. No calcula- 
tion was made for m=0.2 to correspond to the last case 
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under m=0.4 because at the lower Mach number it was not 
possible to calculate satisfactory values of a out to the wing 
tip. The tapered wings were chosen to show, by comparison 
with the first two of the un tapered m=0.4 wings, the effect 
of taper with the span held constant and, by comparison 
with the second and third untapered m=0.4 wings, the 
effect of taper with a given aspect ratio. 

SUMMARY OF COMPUTATIONS 

With m=0A and ^slCc = 0.6j the trailing-edge Mach lines 
do not intersect the leading edge, and the values of Cl^ were 


obtained entirely by means of the conical-flow formulas, as 
follows : 


Component of lift 

Equation 

No. 

Tapered wing, /?.4 = 1 .6 

Untapered, ^.4 = 1.2 

^L!qaC{? 

% total 

^L!qaC(? 

% total 

Uncorrected triangular wing 

Tip effect 

Symmetrical trailing-edge correc- 
tion 

Oblique trailing-edge correction. . . 

(85) 

(96) 

(103) 

(105) 

2.093 

-.190 

-.159 

-.015 

121.1 

-11.0 

-9.2 

-0.9 

2.595 

-.422 

-.340 

-.019 

143.0 

-23.3 

-18.7 

-1.0 

Totals - 

1.729 

100.0 

1.814 

100.0 

fiCi =fiLfqaS, per radian 

1.920 

1.512 


The calculations for tiie remaining values of fiA are summarized in the following table: 


Component of lift 

Obtained from 

Un tapered wings 

Tapered wing 
m=0.4; m»=0.6; 
/te=0.8; /3.4=2.4 

7»=0.2 

m=0.4 

M=0.6 

M=0.8 

PA^l.2 

M=l.6 

fiA=2A 

M=3.2 


Lift on inboard portion 

Lift on outboard portion 

Tip correction 

Fig. 29 or equation (113) 

Equation (115) 

Fig. 31 or equation (118) 

0. 366 
.180 
-.085 

0. 366 
.386 
-.090 

0. 366 
.830 
-.098 

2.128 

.849 

-.363 

2.128 

2.593 

-.392 

2.128 

4.459 

-.415 

1.981 

.462 

-.092 



Totals, ff^LIqac^ 

.461 

.77 

.662 

.83 

1.098 

.92 

2.614 

1.63 

4. 329 
1.80 

6. 172 
1.93 

2.351 

2.20 

^CL^=pLlqaS, per radian 


DISCUSSION OF RESULTS 

The results of the calculations are plotted against the 
reduced aspect ratio fiA in figure 32. The curves for the 
untapered wings may be se en to be approaching, at the 
upper end, the value 2iTirnl^J\ — m^ given by simple sweep 
theory. 

At the lower end, the curves should approach the origin 

along the hne A given by low-aspect-ratio theory 

(reference 13). The two points on the m = 0.2 curve for 
pA<C,l are not entirely accurate because no account was 
taken of the interference between the flow fields from the 
tips. The points are included, however, because, with so 
much sweep, the wing areas affected are small and the inter- 
ference effects should be negligible. The resulting curve 
appears consistent with the corresponding curve calculated 
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Figure 32.— Variation of lift-curve slope with aspect ratio. 


by the slender-wing theory of reference 3, although a dis- 
parity in plan form lessens the significance of the comparison. 

The slender-wing-theory values are also plotted for m=0.4. 
In that case, however, the assumption of extreme slenderness 
is no longer justified and introduces an appreciable error. 
(It should be mentioned that the asympotote for the slender- 
wing-theory curves is below the value given by simple sweep 
theory by the factor V 1 — 

An estimate of the accuracy of the lift formulas of the 
present report, compared, with results which would take into 
account all the successive reflections at the tips and trailing 
edge, may be made from the following observations: 

The values obtained (in the first table) from equations 
(96) and (103), which combine primary and secondary 
corrections, differ from values obtainable for the primary 
corrections alone by only 1 percent of the total lift in the case 
of the tapered wing, and 4 percent of the total lift for the 
untapered wing. Third-order corrections would be only a 
fraction of those small corrections and would, in turn, be 
partly canceled out b}" a fourth-order correction. 

The results in the second table, incorporating the two- 
dimensional approximations, agree within 2 or 3 percent 
with values calculated entirely by the conical-flow method. 

IV— DRAG DUE TO LIFT 

The drag due to lift of a wing with supersonic leading 
edge is simply the lift times the angle of attack. When the 
leading edge is subsonic, the drag is reduced by a suction 
force due to the up wash around the leading edge. In the 
linearized theory this force appears as the limit of the 
product of an infinite velocity across an infinitesimal frontal 
area. 

The formula for the suction force on a subsonic leading 
edge has been derived (see, e. g., Hayes, reference 18) by 
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assuming the flow near the leading edge to be essentially 
two-dimensional and applying the results of two-dimensional 
potential theory. The simple result obtained in that 
manner has been verified for the swept-back wing of finite 
span by application of the somewhat different approach of 
reference 19. 

By the two-dimensional approach, the suction force is 
found to be proportional to the square of the strength of 
the leading-edge singularity in the perturbation velocity u. 
The latter is the quantity discussed earlier in connection 
with the adjustment of the two-dimensional loading to the 
loading on the swept-back wing. With the use of the pre- 
vious terminology it is possible to write for the longitudinal 
component of the suction force per unit stream wise length 
of leading edge, 


dT pTT 

dx m 


V I — 


(119) 


where (equation (53)) is the value, at the leading edge, 
of the coefficient of in Ua. 

Then, if the trailing-edge Mach line does not intersect the 
leading edge, the thrust is merely 






dx 


P TTUq^^V 

2m^ 


Vl — 


( 120 ) 


The total drag due to lift is obtained by subtracting the 
thrust from the product of the lift and the angle of attack, 
or, in coefficient form, 

Cd^ocCj^—Ct ( 121 ) 


where Ct is the thrust coefl5cient TjqS. Thus, in the fore- 
going case, 






( 122 ) 


When a portion of the leading edge is influenced by the 
trailing edge, the leading-edge singularity takes on, for that 
portion, the value given by expression (56), which then 
replaces Ck in equation (119) for the thrust. The total 
thrust is 

J;[a + (AO).+J;’^i«]’*} (123) 

where 


locates the intersection of the trailing-edge Mach line with 
the leading edge. Integrating the first term and reducing 
so coefficient form gives 

so that 


1 Q 

In figure 33, — times the drag-rise factor is plotted 
P Ul 

against the reduced aspect ratio 0A for two combinations of 

sweep and Mach number, m=0.2 and m=0.4, for untapered 

wings. Comparison is made with a theoretical minimum for 

slender wings in supersonic flight obtained by R. T. Jones 

in an unpublished analysis. Using a method similar to Hayes 

(reference 18) and assuming the wing to be narrow compared 

with the Mach cone, Jones has derived a minimum wave-drag 

coefficient 

<‘25) 

where is the aspect ratio defined in the streamwise, 
instead of the spanwise, direction; that is, if I (numerically 
equal to Xt) is the over-all length of the wing in the stream 
direction, 

Ar=PIS (127) 
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Ficure 33. — Variation of drag-rise factor with aspect ratio for untapered wings. 


The wave drag is to be added to the vortex drag, which is 
the induced drag of subsonic flow, calculated from the 
spanwise loading. Using the minimum induced drag ob- 
tained from lifting-line theory gives as the minimum super- 
sonic drag-rise factor ® 


1 , 

Cl TrA~^2TrAx 


(128) 


It may be seen that the drag rise of the constant-chord 
swept-back wings is fairly close to this minimum, especially 
at the lower values of m for which equation (128) was 
derived. 


V— SUMMARY OF FORMULAS 

The formulas for the loading, lift, and drag coefficients are 
summarized in the following table, in which the equations 
are identified by number. 


* This result has sinoe been published in The Journal of the Aeronautical Sciences, vol. 
18, no. 2, Feb. 1961, pp. 76-81. 
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Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Mar. 15, 1950, 
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APPENDIX A 

SYMBOLS 


V 

M 

fi 

U, V, w 

p 

5 

Ap 

a 

L 

T 

Cl 

Cd 

Ct 


Co 

Ct 

s 

S 

I 

A 

X 

A 

A. 




Xa, Va 


XbyVb 

XcjS 

XofVo 

XuVi 

X2yV2 

x*,y* 

XtyS 


GENERAL 


free-stream velocity 
free-stream Mach number 


VM^-1 

perturbation velocities in streamwise, cross-stream, 
and vertical directions, respectively 
density of air 


dynamic pressure 



pressure difference between upper and lower 
surfaces, or local lift 
angle of attack, radians 
lift 


leading-edge thrust, or component of leading-edge 
suction force. in flight direction 

lift coefficient 


(A) 

lift-curve slope 
drag coefficient 
thrust coefficient 


WING DIMENSIONS 

root chord 
tip chord 
semispan 
wing area 

over-all length in the steamwise direction 
angle of sweep of the leading edge 
taper ratio (cjco) 
aspect ratio (4iS^/S) 
streamwise aspect ratio (PIS) 


RECTANGULAR COORDINATES 

Cartesian coordinates in the stream direction and 
across the stream, in the plane of the wing 
coordinates of apex of conical field used to cancel 
triangular-wing loading (Equation (8) at tip, 
equations (21) and (22) at trailing edge) 
coordinates of apex of conical field used in second- 
ary cancellations 

coordinates of point on tip; apex of conical field 
used to cancel assumed cylindrical load 
coordinates of intersection of Mach forecone from 
x,2/ with edge at which correction is being made 
coordinates of intersection of trailing-edge Mach 
cone with leading edge (xt given by equation 
(61)) 

coordinates of intersection of Mach line from Xi, yi 
with trailing edge (x 2 given by equation (112)) 
coordinates of intersection of tip Mach line with 
trailing edge 

coordinates of intersection of tip and trailing edge 


? streamwise distance of x^y back from leading edge, 

as a fraction of the tip chord (equation (69)) 
fo distance of Xq,s behind leading-edge tip, as a frac- 

tion of the tip chord (equation (70)) 
fc distance of XcyS behind leading-edge tip, as a frac- 

tion of the tip chord (equation (64)) 


CONICAL COORDINATES 


ao 

CL 2 

di 


To 


In the following, all slopes are measured counterclockwise 
from a line extending downstream from the apex of the wing 
or of the pertinent canceling sector: 
slope of leading edge 
slope of Mach lines 
slope of trailing edge 
slope of Mach lines 
slope of ray from the origin _ y 
slope of Mach lines ^ x 

the value of a corresponding to a primary canceling 
element of which the apex lies on the Mach fore- 
cone of the point at which the load is being 
calculated (equation (13) for tip corrections, 
equation (25) for trailing-edge corrections) 
limiting value of a for leading-edge correction 
(equation (47)) 

^( 2 ^ 2 , 2 / 2 ) (equation (108)) 
a(x„ s) 

slope of ray from apex of element d_ay—ya 


m 


m 


do 


-=)3 


X — Xa 




y—s 


slope of Mach lines 
slope of ray from x„,y„ y—y^ 
slope of Mach lines ^z—Xn 
slope of ray from XcfS__ y^-s 
slope of Mach lines ^ x — Xc 

slope of ray from leading-edge tip 
slope of Mach fines 

slope of ray from ^y — y* 

slope of Mach fines ~^x — x* 
limiting value of ^0 for leading-edge correction 
(equation (39)) 

limiting value of ta for leading-edge correction 
(equation (44)) 


x — {Pslm) 


COMPONENTS OF STREAMWISE PERTURBATION VELOCITY 


u^ basic (uncorrected) perturbation velocity as given 

by solution for triangular wing (equation (6) 
for subsonic leading edge) 

Uq value of at a=0 (equation (7)) 

Au correction to induced by cancellation of pressure 

differences outside the wing plan form 
Ua constant perturbation velocity on sector used in 

canceling triangular-wing loading 
Ut, constant perturbation velocity on sector used in 

secondary cancellation 

Uc constant perturbation velocity on sector used out- 

board of tip in canceling assumed cylindrical 
field 
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(Ai^)o 

(A2U)o 

(^U)a 

Au* 


symmetrical trailing-edge correction to Ua (equa- 
tion (20)) 

correction induced by canceling (Au)o at leading 
edge (equation (38)) 

correction to Ua due to single oblique trailing- 
edge element (equation (24)) 
value of tip (’orrection to Ua at the point 


ARBITRARr MATHEMATICAL SYMBOLS 


(AOo 

dAC 


da 


da 




M 

9 

C 

R 

r. p 


value of coefficient of , ~ = ^ in Ua at the lead- 
^mx—rpy 

ing edge (equation (53)) 

decrement in. Ca due to reflection of {Au)o at lead- 
ing edge (liquation (54)) 


decrement in Ca due to reflection of (Au)« at lead- 
ing edge (equation (55)) 

Non-dimensional expression for strength of the 
leading-edge singularity (equation (59) ) 


value of <T at leading-edge tip 



taper parameter I — ^ — j 


function defined by equation (79) 
inverse-cosine term of leading-edge correction 
function (equation (35)) 

radical term of leading-edge correction function 
(equation (36)) 
real part 


ELLIPTIC INTEGRALS AND FUNCTIONS 

k modulus of elliptic integral, defined where used 

(also with subscripts) 
complimentar}" modulus ( 

<t> ov \l/ argument of elliptic integrals, defined where used 
(also with subscripts) 

F((t>,k) incomplete elliptic integral of the first kind of 
modulus k and argument <t> 

K,K(k) complete elliptic integral of the first kind; that is, 

t) 

E(<t>jk) incomplete elliptic integral of the second kind of 
modulus k and argument </> 

£*, E{k) complete elliptic integral of the second kind ; that is, 

E=E(^.i) 

Ko -K 

IT 

Eo -E 

7T 

K' K{¥) 

E' E{k') 

Z zeta function (equation (41)) 

Ao function used in evaluation of elliptic integral of 

the third kind, circular case (equation (16)) 

fi function used in evaluation of elliptic integral of 

the third kind, circular case (equation (Bll)) 


APPENDIX B 

EVALUATION OF THE INTEGRAL IN EQUATION (26) 


It is first necessary to recall that ta is a function (equation 
(23)) of X, y, and a. After substitution for (« in equation 
(26), we may integrate by parts to obtain 


f “0 

Jo 


_1 (1 — o) (Ig— wi,)— (w<— g) (1—0 du^ (g) 


— «o 


(1— m,) (<a— a) 

_j (1 — 2w,(x— Co) 


da 


cos“‘ — 




(x 


-Py)-yj 


fiy—m,(x—c^ 

x—^y—m,Co 


X 


r 


da 


iPy—ax) V(l— a) {ao—a) {m—a) (m+a)J 


(Bl) 


The integral term on the right-hand side of equation (Bl) 
is an elliptic .integral of the third kind which may be eval- 
uated through the substitutiou of 


U3=sn 


2 m (ao— g) 
(m + ao) {m—a) 


, __ / (I— m) {m + ap) 
2m(l— ao) 


If the value of co at the lower limit is designated by wo, this 
substitution gives 


Jo (By- 


da 


(Py-ax) ^{l—a) {ao-a) {m—a) (m+al 




A CXZ f 

— aoX V m(l— ao) Jo 


, m + ao_2 

, 1 

"0 2m 


where 


n= 


or 


rao 
Jo I 


1+71 sn^o) 

( m + gp) {mx — fiy) 
2m{fiy—aox) 

da 


(B2) 

(B3) 


0 {py — ax) ^J{l—a) {ao—a) {m — a) (m + a) 


mx 

where 


, , . f"0 d03 

n3(o)o,A:,n)-J^ 1-1-M.sm'a> 
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is the normal form of the elliptic integral of the third kind. 

It is first noted that n>0. For this case it can be shown 
that the substitution 


f) 


(B5) 


gives 


X. / IN 2/ 7 A I sn{v,k')cn{v,k')r wo . , ^"1 

— uyoCU {vjk)'{- du(y k') 

(B6) 


where 




(B7) 


is the amplitude of the elliptic function v, Ao is the function 
defined in equation (16), and 12 is an angular function of k,v 
and too which will be discussed later. 

If 


, • -1 / 2tto 

i=sm *-»/ — |— 

V m + Oo 


then 


u>o=F{>l/,k) 


From equation (B5), sn{v, k')=-^ n+P ’ k ') — 


(B8) 

(B9) 

~P~ 




Vk^d +n) 

may be found, SO that equation (B6) 

may be rewritten without recourse to the Jacobian elliptic 
functions as 

IlaCtoo, k, Fi'f', k) + 


^ {l+n){n+k^\_ ^Ko 

This expression is to be substituted in equation (B4) and 
the result used in equation (Bl). As previously mentioned, 
the functions Ko and Aq are tabulated in reference 1 1 and 
is plotted in figure 6. The function U is given by ^ 


22 (- sin 2j^ sinh 2j~ 

i=l fis ik. 


12= tan 


1 —222 (— cos 2j ^ cosh 2j^ 


(BID 


with 


q=e K 


(tabulated in reference 15). 


® The S 3 fmbol q in equation (Bll) is standard notation for the nome of the Jacobian theta 
function, and is not related to the dynamic pressure q of the text. 


APPENDIX C 

INTEGRATION FOR LOSS OF LIFT AT THE TIP OF WING WITH SUBSONIC LEADING EDGE 


From equations (88a) and (6) 

< UL\ f"* G'{ a) 

J«< 


(-) - 

\q.oc/tip 


da 


(Cl) 


where 


Va Jat 

C>( \ [{at . m,—a, \/ I o+F a \ 

a,*a*(TO<— a)L\® m,— a/ \V + m,/ 

2 Vw,+m<*V®+® J ^ 

The terms in G'{a) are of two types; namely, those that 
contain -^a+a^ and those that do not. The fornjer combine 
with the radical ■yjrn^—a^ in equation (Cl) to form elliptic 
integrals of the first, second, and third kinds. The latter 
give rise to terms in equation (Cl) which are integrable by 
elementary means. It is convenient, therefore, to consider 
the integral in two parts, writing 


_ r 


^Jm' 

where h is that part of the integral not requiring elliptic 
integrals. 

Then 


a— a, 1 

/at . m, — aA 

, da 

2m(l+m) dll 3 (w,,fc,n)“|) 

jat aYmia{mt — (iy 

\a rrit — aj 


m,—m bn Jj 


1 


r {mt — atY 




_i m<at— 


mtat — m^ 




-af^ 


(C3) 


The remaining terms, involving + ^nd 
are integrated by means of the substitution 


The result is 
/2= 


sno) 


1 


V m — a j __ j\ 

m(l + a) 


— m 


(C4) 


aYmt^|2mmt{\-\-m 
at{m — at) 


(r mt — at / rrit — at at\ 
^(L rut \l+mt m) 






m,— a,\ Im + a, . 

(“ — sr) ^ 


(l+wt)(m.— a.)^ 




( 05 ) 
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where cvt = w(at)y \}/ is its amplitude, and 


n=~ ^ 

rrit — m 

(C6) 

From equation C(4), 


. • >1 / m — at 

i/'=sin — T 

^7n{l+at) 

(C7) 

The elliptic integral 


Ui{u>„k,n)= f“‘ , 

^ ^ Jo 1+nsn^co 



is evaluated in equation (BlO). Its derivative with respect 
to the parameter n may be obtained for this case (n^O) in 
the form 


51I3 

dn 


2n ^l 0 


n 


2n («+F)(l + 71 ) 

^)J— -y/j 




S+F)(TTS) L^^Cp+I+p) 


where 0 and are the angles defined in equations (B7) and 
(Bll) and the elliptic functions cncvt and dncot, obtained 
from equation (C4), have the values 


cno)t = 


V at{l+m) 
m{l+a,) 


dno)t= 


V (1 -|-w)(m-|-a() 
2m(l 4-a,) 


(C9) 
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